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EXECUTIVE SUMMARY 


The ongoing effort of introducing wireless communication systems onboard naval 
assets provides a noticeable incentive for developing software tools capable of predicting 
the signal strength variations in indoor environments. To this end, this thesis focused on a 
tool developed by the Ohio State University Electroscience Lab (OSU-ESL) called 
Numerical Electromagnetic Code-Basic Scattering Code (NEC-BSC) and evaluated its 
applicability to confined spaces. Initially, some simple geometries were created in an 
effort to gain familiarity with the code and its capabilities and limitations. Even at this 
early stage, it was found that the code erroneously generated electromagnetic waves 
penetrating the surface of a perfectly electrical conducting (PEC) cylinder. However, a 
comparison of NEC-BSC predicted values with the results presented by a different 
approach called shooting-and-bouncing ray method (SBR) in a relatively simple 


rectangular room with furniture offered encouraging results. 


The major task, in any event, was to compare the signal strength predicted values 
obtained by NEC-BSC to actual measurements conducted inside a building at the Naval 
Postgraduate School in the 2.4 GHz ISM band. The comparison was executed in two 
ways. First, by contrasting the NEC-BSC and physical electric field strength coverage 
maps and, second, by creating two data matrices based on the predicted and measured 
values. In the former case the two maps agreed very well, whereas in the latter the 
agreement was not illuminating and that could possibly be attributed to the many inherent 
errors associated with the second approach. In parallel, a statistical analysis of the electric 
field strength variation was performed in the same building using NEC-BSC. In 


particular, the antenna was placed in five different locations and, subsequently, five sets 


AV 


of data were acquired. Then the excess electric field “room gain” over free space was 
calculated for all same-floor points within certain distances from the antenna called 
“distance-groups. In this manner, a random variable was defined describing the variation 
of the so-called “room gain” inside the building under consideration and its probability 
density function (pdf) and cumulative distribution function (cdf) were computed. It was 
found that regardless of the distance between the transmitter and the receiver on the same 
floor the aforementioned random variable had a lognormal distribution and its pdf and 
cdf obeyed a Gaussian pdf and cdf, respectively, having identical mean and variance. 
Furthermore, the modeled floor of the building was found to exhibit gain over free space 
for all antenna locations. Finally, the mean value of the “room gain” random variable 
ranged from 4.276 to 8.322 dB while its standard deviation varied from 4.624 to 5.275 


dB. 


XVI 


I. INTRODUCTION 


This highly technological era has been called “digital revolution,” in 
; I ae - , 
correspondence with the 19" century's “industrial revolution.” An important long-lasting 


contributor to this revolution 1s the wireless Communications. 


A. MOVING TOWARDS WIRELESS 


Wired communications have offered for quite a long time a reliable and speedy 
means of exchanging ideas and promoting cooperation amongst businesses and 
individuals. However, as the demand for more capacity and flexibility grew, new 
solutions had to be considered. “Mobility” was the new logo and that entailed terminals 
and equipment moving around, installation being easy and quick, inherent resistance to 
natural and manmude disasters being vital. The cost of uchieving all that with 
conventional approuches was posing an impediment to progress. Wireless 
communications not only satisfied the above requirements but they also introduced new 
ideas like roaming, creative solutions like setting a network at locations where cable 
could not go, and modern schemes which could support already existing technologies like 
the Internet. Furthermore, wireless communications helped in cutting down expenses 
since the exorbitant price of replacing the wiring in and around a building was no longer 
an issue, promoted the concept of cells (e.g. megacells, macrocells, microcells) thereby 
enubling communications at high rates, and made communications available to larger 


poruons of the population. 


Bb. INTEREST FOR WIRELESS INDOOR USE 


Nevertheless, there were still persons and sites excluded and a huge increase in 
the number of potential users. Customers were experiencing a shift to yet another idea: 
service provided anywhere. anytime. Moreover. the planning for broadcasting services 
hke Terrestrial Digital Audio Broadcasting (T-DAB) or Digital Video Broadcast 
Television (DVB-T) [1] conveyed a clear message: Picocells had to be invented and 
coverage in confined environments should be taken for granted. Additional motivation 
was given with the advent of Personal Communication Systems (P.C.S.), and various 
emerging wireless standards like IEEE STD 802.11 [2], HIPERLAN/2 (High 
Performance Radio Local Area Network) [3], HomeRF [4], and Bluetooth [5]. 

Until recently, the majonty of indoor wireless systems occupied the 1.8-2.0 GHz 
band (e.g. Digital European Cordless Telephone-DECT, Personal Access Communication 
Systems-PACS. Personal Handyphone System-PHS) {6]. Yet, this trend has shifted to the 
license-tree Industrial Scientific & Medical (ISM) bands, 902 to 908 MHz, 2400 to 
2483.5 MHz, and 5725 to 5850 MHz. These bands allow fee-free operation, inasmuch as 
no one actually owns these frequencies. The 2.4 to 2.485 GHz band attracted some extra 
preference on the part of the business world. Europe and the US adopted the same 2.4 
GHz band. whereas Japan allocated only the frequencies 2.471 to 2.497 GHz. In addition, 
apart from some excepuons like Nokia Rooftop [7] operating at 5.8 GHz, most of 
companies have largely invested in 2.4 GHz projects. The reason appears to be two-fold. 


First. the bandwidth of ISM 2.4 GHz is three times that of 900 MHz, and, second, 2.4 


fh 


GHz 1s subjected to less “radio-traffic” and its resulting interference. There is also move 
into the 5.2 GHz ISM band because of un even wider bandwidth. 

Remarkable and tempting as all the aforementioned arguments may be, they fail to 
reveal the difficulties of indoor propagation. Simply put, radio transmission and reception 
within buildings or other enclosed environments suffers from various phenomena. 
Primarily. Multipath fading 1s caused by the simultaneous arrival of electromagnetic 
(EM) waves from different directions via reflection, refraction, diffraction, and scattering, 
which add either constructively or destructively. Consequently, the received signal varies 
dramatically. even if we move the mobile antenna as little as a half-wavelength distance 
In any direction. In this context, the signal strength may degrade severely the 
communications range. 

Despite the problems associated with the indoor radio channel, the necessity for 
meeting the demand for microcellular and picocellular communications remained. 
Research on this field was launched as early as 1982 by Bnitish Telecom [6]. 
Meusurements at vanous frequencies inside diverse buildings were performed. Soon, a 
new requirement for the determination of the link budget was realized. Should a robust 
Wireless System be designed, the modeling of the indoor radio channel had to be done 
successtully. 

Essentially. until the mid-90’s two kinds of indoor propagation models were being 
ulllized: Empirical and Physical. The former 1s bused on actual measurements run in the 
locations of interest with the intention of coming up with an analytical formula, which 


could later be applied to a similar environment. The latter depends on the laws of physics, 


the UTD/GTD (Uniform and Geometrical Theory of Diffraction) in particular, and uses 
site specific (SISP) information. Both models strive to determine the coverage area, 
transmitted power requirements, number of base stations or access points (AP) (Figure 
1.1), optmum locations for antenna mounting, and all other major concerns of 
communication engineers dealing with indoor wireless communications system (W.C:S.). 
Further discussion on indoor radio propagation and the models pertaining to it will be 


made in the following chapter. 





Figure 1.1: Access Point Installed at the Network Lab of the NPS ECE 


Department (From [29]) 


oa DEVELOPING A SHIPBOARD W.CSS. 


Considering the numerous advantages of W.C.S., the U.S. Navy recently launched 
a thorough study on this field [8,9,10,11]. In 1999, test, evaluation, and installation of 


personal communication system (P.C.S.) onboard various U.S. ships took place [12]. 


Although research 1s still in progress and the market has not arrived at a mature Stage, the 


benefits of implementing a shipboard W.C.S. are pronounced and well established: 


Availability: Current commercial solutions for wireless local area network 
(WLAN) may already satisfy the Navy’s requirements [11]. Selection could be 
made over a broad variety of packages so long as they are compatible with the 
[BBEest® S021. 

Affordability: Not only the acquisition cost will be low but, by virtue of cabling 
elimination, the support and maintenance expenses will be reduced as well. The 
W.C.S. could be built with ruggedized commercial off-the-shelf (COTS) digital 
equipment rather than Navy-proprietary solutions. 

Flexibility and Mobility: The W.C.S. could be easily installed in previously 
unreachable locations on the ship. Furthermore, crew members could have access 
to a centralized system from virtually any ship compartment and, in tum, their 
whereabouts could be monitored providing vital information, should any 
dangerous situation arise. 

Survivability: Wireless systems mitigate the effects caused by cable plant 
damage, fire or other minor or major disturbance. A fault-tolerant operating 
system should be the ultimate goal. 

Back-up: When mostly needed. in case of emergency. WLAN could either be 
used for damage control communications or simply substitute the existing wired 
LAN. For this purpose, an anti-fading system with a directional antenna could be 


implemented to account for fire and the fire extinguishing effects [13]. 


e Improved data handling: Information collection will be done expeditiously and 
report assessment will be carried out on the spot, enabling instant digital 
processing. 

e Expandability and variability: Having already been tested for interoperability in 
the commercial sector with encouraging results, WLAN’s are expected to be 
Incorporated in various arrangements and offer an architecture which will support 


a plethora of configurations without sacnficing the end-to-end performance. 


D. GOAL FOR THIS THESIS 


Driven by the ongoing effort of introducing W.C.S. onboard Naval assets, this 
thesis addressed the problem of indoor propagation modeling using a tool] developed by 
the Ohio State University Electroscience Lab (OSU-ESL); specifically, the Numencal 
Electromagnetic Code- Basic Scattering Code (NEC-BSC). The approach initially 
consisted of creating simple geometnes in order to gain familianty with the code and its 
capabilities. The next was the simulation of the second floor of Bullard Hall, at the Naval 
Postgraduate School (NPS), in Monterey, California. Finally, a comparison was made 


between the NEC-BSC predicted results and those measured M.M. Matthews [11]. 


E. THESIS OUTLINE 


This thesis is divided into seven chapters. The second chapter discusses the 
theoretical background for an indoor radio channel and gives insight into previous work 
on predicuon of radio propagation in confined environments. Chapter [IJ introduces the 


NEC-BSC and summarizes its limitations and capabilities. Next, vanous simple 


simulations involving uncomplicated geometries are run and explained in Chapter IV. 
Chapter V describes the modeling of the simulated floor and the simplifications adopted 
in our model. All results and the statistics extracted from the simulations are shown in 
Chapter VI. Finally, in Chapter VII, the conclusions are presented and some 


recommendations for further work are made. 
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I]. BACKGROUND 


A. THE INDOOR RADIO CHANNEL 
by Complex Environment 


The indoor radio channel 1s a complex propagation environment and there 1s good 
evidence through years of research that it 1s more complicated than an outdoor radio 
channel. This is due to the vast diversity of the potential surroundings. 

Broadening the service area to reach more end users at higher data rates and 
virtually everywhere was not an easy task to undertake. Eventually, picocells had to be 
designed and put into practice. In doing so, businesses extended wireless networks 
coverage to include office buildings, industrial plants, railway stations, tunnels, and 
airports. Users could be either static (e.g. by establishing a WLAN between PC’s) or 
mobile (e.g. by expanding cellular telephony services). 

Prior experience obtained in various outdoor schemes does not guarantee 
applicability in enclosed sites. Unlike city and rural topologies, the indoor environment 
consists solely of man-made structures and obstacles like walls, partitions, windows, 
furnishings etc. In addition, in indoor cells there 1s no material uniformity, and many 
more flat surfaces compared to outdoor cells. Given that the frequencies of interest reside 
in the ISM band, the dimensions of obstacles are typically much larger that the radio 
wavelength (2). Moreover, the observed channel depends on the type of building, whether 
there are doors opened or closed, location of antennas, presence of people. etc. 


Shadowing regions are also an issue because by placing the transmitter and receiver 


inside the same buijding, many obstructed paths are created and diffraction tends to be a 
significant contributor to the received signal strength. Even the existence of corridors 
plays a role since it exhibits a waveguide behavior [14,15]. All these parameters have a 
nouceable effect on the resulting propagation. They produce a multipath environment in 
the sense that many potential transmission and arrival paths can exist. 

As opposed to macrocells, picocells, sometimes called indoorcells, have been 
found to be more deterministic than statistical and are not stationary in space and time 
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Figure 2.1: Example of Multipath Fading line-of-sight Case (From Ref [17]) 
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Figure 2.2: Example of Multipath Fading Obstructed Case (From Ref [17]) 


ize Multipath Fading 

In this complex environment, the EM waves generated from the transmitter reach 
the receiver Via numerous paths with unequal contribution. A line-of-sight (LOS) may not 
exist (figure 2.1), in which case only obstructed (OBS) paths are present (figure 2.2). The 
signal traveling via these paths can be gravely attenuated because of wall penetration 
and/or reflection/diffraction from various absorbing materials on the premises. However, 
even if a LOS ts present there are stil] other propagation paths whose contnbution might 
be of lesser significance yet not at all negligible. Consequently, the signal arnving at the 
receiver will be distorted with the sum of all signals received tending to have random 
phases. The received signal power usually vanes as much as 30-40 dB (small-scale 
fading) when the receiver 1s moved for only a fraction of wavelength. Statistical analysis 
has shown that the magnitude of the sum of arriving signals may be described as Rician 
random variable (when a LOS path exists) or as Rayleigh random varialbe (if only OBS 
propagauon paths are present) [17]. In particular, the phenomena leading to multipath 
fading are the following: 


sea 


a) Reflection 


Reflection comes about when an EM wave strikes an obstacle with 
dimensions much larger than the wavelength. In this case, the surface is said to be 
“smooth” and the wave specularly reflected. The criterion used in deciding whether this 


behavior should be anticipated ts the Ravleigh Criterion: 


Lf oa (2218 
8cos @, 


where Ah is the height difference between two points on the surface under 
consideration, 4 the wavelength, 9 the incidence angle. If the inequality (2.1) holds, then 


the Rayleigh criterion is met and the surface 1s assumed to be “smooth.” 


b) Scattering 


Scattering happens when the object hindering the EM wave's path has 
dimensions either comparable to or smaller than the wavelength. The EM energy is 


reradiated in all directions as shown in Figure 2.3. 
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Figure 2.3: Effect of Surface Roughness on Reflection (From Ref [17]) 


C) Refraction 
Retraction represents the change in the direction of propagation of an EM 


Wave as It crosses the interface of two media with different velocities of propagation. 


d) Diffraction 
Diffraction occurs when an obstruction made of material impervious to 
EM energy stands in the path of propagation. The area behind this obstruction is defined 
as the Shadowing Region. Huygen’s principle in this case dictates that secondary EM 
waves will be generated from the edges of the obstacle and EM energy will eventually 
“creep around’ and propagate in the shadowing region, albeit with reduced magnitude. 
Thus, diffraction contributes substantially to signal reception in non-line-of-sight (NLOS) 
situations. 
a Path Loss 
In an effort to design and establish a W.C.S. in an indoor environment, the 
communications engineer will need prior accurate knowledge of certain factors like the 
potential coverage area. best locations for antenna placement, transmitting power 
requirements, etc. In short, the link budger, in which all gains and attenuations of the 
channel are considered with the ultimate goal of estimating the performance of the 
communication system in the given surroundings, needs to be determined. An essential 
tool indicative of the expected communications link quality 1s the path loss (PL), which 1s 
the difference. expressed in dB, between received and transmitted power. An example of 


what can be expected in an indoor office environment is given in [18], where 


i 


measurements at 2.4 GHz were performed and it was found that the maximum acceptable 
PL was 123 dB with the SNR (signal-to-noise ratio) always more than 20 dB. 

Early works in this area revealed an inclination to a simplified approach in 
describing the indoor channel. It was thought that the PL models used in the mobile 
channel could just as easily be conveyed to indoor applications by adding an attenuation 
factor (AF) to account for walls. floors, etc. The attenuation factor model was a paradigm 


In this regard: 


PL Gigs | — Pita ids alone los) + FAF[dB] (2.2) from [6] 


0 
where d and d, represent the distance in meters from the transmitter to the 
receiver and to a reference point. respectively. 2,, denotes the attenuation exponent for 


the “same floor” measurements, and the FAF is simply a factor whose values usually 
depend on the radio frequency and the number of penetrated floors. 
In the same fashion, the Keenan model takes into account the number of floors 


and walls between transmitter and receiver, nm, and 1,,, respectively, and the 
corresponding attenuation factors, 1.e. a, and a,,: 
b= 1 20Rse 4 1t) Greriat, ay (2.3) from [17] 
where 1, denotes the loss at distanec = Im: 


However, this concept 1s not justified in [14] where it is argued that the resulting 


error margins reach unacceptable levels in some indoor environments. 


Indeed [19] and [20] clearly point out that the PL within buildings does not 
depend on the distance alone, but also on the floor area, the number of walls, partitions 
und obstacles between the transmitter and receiver, and on the signals leaving and 
returning to the building (in case of proximity to neighboring buildings). In addition, the 
survey in [14] concluded that the frequency of the radio wave is also significant inasmuch 
as the higher it 1s the lower the penetration, and thus the path loss. Through a different 
perspective, in [2]] it 1s suggested that another concern is the path the EM wave follows 
to reach the receiver: specular paths result in cumulative transmission PL whereas the 
diffracted ones (whenever not passing through many walls) display a lower PL. 

It should be stressed, nonetheless, that there 1s considerable difference in the 
observed staustics of PL between the “same room” and “many rooms” propagation case. 
In the former case the PL exhibits a slow-fading fluctuation profile, but at the same ume 
the attenuation factor approaches the free-space n = 2 [22]. In the latter case, on the other 
hand, the further away the transmitters room is, the attenuation factor may take very 
different values. For instance. in [23] and [24] it was estimated that n = 6.0 at 2.3 GHz 
and in [14] it was estimated that 7 = 3.0 for an office building. 

A typical example of the immense diversity characterizing the results reported in 
the literature is given in [14] where the measured PL through concrete ranged from 7 dB 
to 27 dB depending on the approach chosen. By contrast, in [15] it was estimated at 3.5 
dB. 

Although it is true that PL depends on the number of obstructions (which are 


many in the indoor environment) between transmitter and receiver, it 1S also true that not 


all of them are deemed significant. An interesting theoretical rule of thumb is whether 


the obstruction in question occupies 0.6 times the first Fresnel zone. In this case, the 


Impact on the received signal strength is low. The radius r, of the n” Fresnel zone is 


given by (assuming r, << d,,d,): 


ddd, 
S| aeeee (2.4) from [17] 
ed. 


All the pertinent distances are shown in Figures 2.4 and 2.5. 
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Figure 2.5: 0.6 times First Fresnel Zone (From [17]) 
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B. MODELLING THE INDOOR RADIO CHANNEL 


Given the increasing demand for indoor W.C.S., installation in diverse confined 
environments and the driving force of business competition dictates a prompt and reliable 
mechanism of predicting the performance of any system under such harsh operating 
conditions. Measuring the signal strength in one building and subsequently estimating 
the best antenna locations appears to be a tempting and simplified approach. However, 
this method would have applicability only in identical buildings since it is well known 
that only slight differences in the geometry of the building or in the properties of the 
materials constituting the surroundings, or even a small alteration in the transmitted 
frequency, would create a very dissimilar situation. Furthermore, performing 
measurements in all existing buildings or other constructions (factories, train stations etc) 
is un unrealistic method, inasmuch as it would necessitate vast human and equipment 
resources. The answer was given by the indoor radio channe] models. 


1. Categories of models 


Extensive research in the literature has indicated that we can distinguish three 


culegories of models: 


a) Empirical 

These are also called statistical models and their principal feature 1s the 
extracuon of an expression, which could be a function of several parameters (e.g. 
ullenuauion factors for concrete walls. windows etc), accurate enough to fit curves of a set 
of measured datu. This expression usually has a broad field of applicability, but cannot be 


universally used nonetheless. Often this means that these models yield satisfactory results 


for a group of buildings that display characteristics similar to the ones where the 
meusurements were taken. Examples of empirical solutions are the Keenan factor model, 
the COST231 {17]. the Encsson Breakpoint model [6], and the model of [18] where 


attenuation due to walls was not taken into account. 


b) Physical 

These are also called deterministic or site specific (SISP) models. They 
have recently dominated the field of indoor radio prediction tools since their evolution 1s 
dependent on computer growth and improvement. They basically utilize various ray- 
tracing techniques and largely depend on the EM laws of physics and the data of the 
specific site where they are applied. Physics plays a significant role through the use of the 
theory of UTD/GTD and the Method of Moments (MoM) [25]. Moreover, the data 
frequently include the constituuive parameters of the materials found in the site under 
consideration, the detailed topography obtained by blueprints, and the precise location of 
obstacles. A different but useful approach is proposed by finite difference time domain 
(FDTD) models where fewer resource requirements have to be met compared to 
UTD/MoM models [26]. 

SISP models mostly take advantage of the high frequency approach, which 
enables the simulation of radio propagation in the form of rays. Therefore, the resulting 
total field strength at the receiver's position 1s computed as the sum of all the arnving 
rays generated by the transmitter. The rays follow multiple paths consisting of reflections, 
refractions. diffractions, scattering, and also double or triple bounce components (e.g. 
reflecuon/diffraction. or reflecuon/diffraction/retlection). However, not all SISP models 


18 


account for all possible combinations. More specifically, diffracted rays were not 
examined in [19], [21], and [27] while scattering was not considered in [28], [27], and 
[1]. Furthermore. the FDTD model of [26] does not include ground reflections and the 
effect of antenna height. On the other hand, PlaceBase tool takes into account all the 
above but fails to consider furniture and obstacles [29]. 

Examples of GYD/UTD and MoM models are described in [28], [1], [30], 
and {31!] whereas examples of various ray-tracing techniques are analyzed in [19], [32], 
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Cc) Other Models 


There are models that do not fall into the above categories. Hybrid models 
occupy a middle ground between empirical and deterministic approaches. They combine 
the simplicity of the former with the flexibility of the latter producing somewhat efficient 
results. Two best known hybrid models are found in [1] and [35]. An entirely different 
solution is suggested by Artificial Neural Networks (ANN) models, which utilize the 
principles of Feedforward Neural Networks. An example is given in [36]. 


me Language and Input Data Tools 


Currently. as au survey of the literature has indicated, there are at least three 
languages used for code development: FORTRAN (e.g. for NEC-BSC and [30]), C++ 
(e.g. for [33] and [37]), and the Eclipse constraint-logic programming language [34]. 

With the advent of effective software tools. the time and effort for inputting data 
in the model is substantially minimized. Indeed, this is an important factor whenever 


different codes are weighed against each other since not all of them incorporate this 
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capability. For instance, the models of [33] and [35] make use of AutoCAD in order to 
convert the building data into ASCII files. By virtue of these files, the layout of the indoor 
environment is converted into numerical coordinates subsequently used for the simulation 
culculations. By contrast. the alternative to this convenience is to manually input all 
coordinates which. needless to say. 1s a relatively cumbersome task (e.g. NEC-BSC). 
However, we should underscore the fact that AutoCAD files are not readily available for 
all existing buildings. It would be rather optimistic for one to expect to find AutoCAD 
files for a building built in 1950's. 

a Extracted Data 

The extracted results may be given in various ways. The majority of models 
generate ASCII raw files which include important factors such as impulse response, path 
loss. suggested location for AP, or all of them. In particular, PL can be expressed as a 
ftuncuon of either frequency or distance between transmitter and receiver. Again, not all 
models convert automatically these raw files into graphs and the user’s intervention 1s 
required in such cases. Also, some applications, like for instance WiSE, CINDOOR, 
WinProp, SitePlanner, PlaceBase (figure 2.6), initiate a colored graphical representation 


of the overall radio propagation for the building under consideration [29]. 
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Figure 2.6: Bullard’s 2° Floor Relative Signal Strength Coverage Map Generated 


by PlaceBase (From [29]) 


4. Comparison of Empirical & Physical Models 


In what follows, we attempt to contrast the two more popular models by showing 


their advantages and disadvantages. 


ih) 
[4 


a) 


b) 


Advantages 


(1) 


(2) 


Empirical models 


Simple and fast 

No requirement for detailed determination of the 
confined environment. Only basic facets are needed 

All environmental radio propagation factors are 
implicitly taken into account 

Computer specifications for implementing empirical 


models are on reasonable levels 


Physical models 


The more SISP data we include in the model the better 


the predictability for the indoor radio channel 

Reveal more physical details 

[f employed properly, they are more accurate 

Predicted results are achieved without performing 
measurements 

Since they are based on the principles of physics, they 


can be used in a plethora of enclosed sites 


Disadvantages 


Ce 


Empirical models 
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Do not take advantage of possible prior knowledge of 
the physical environment 

Do not provide reliable results if applied on 
complicated indoor situations (e.g. many obstacles, 
people etc) 

Do not take into account the waveguide effect that is 
profoundly important in corndors and tunnels 

Are inextricably dependent upon the accuracy of the 
measurements. All possible errors will affect the 


reliability of their predicted data. 


Physical models 


Any kind of distance or wall characteristics tolerances 
will have a negative impact on the results. 

In case they employ GTD/UTD methods, they 
inherently assume that the surfaces of the buildings and 
the ground are perfect flat planes, thereby neglecting 
scattering caused by small objects attached to those 
Surfaces. 

In general, they are more time-consuming 

They largely depend on the accuracy of the input 


material properties 


ZS 


e They do not yield reliable results for distances far away 
from the transmitter due to the large number of 
interactions required in such cases. Hence the hmitation 
of their coverage area. 

e Since the walls and floors are inhomogeneous, some 
sort of trade-off becomes imperative. As a result, errors 
are likely. 


>; In Which Group Does NEC-BSC Fall? 


Clearly. since NEC-BSC obeys the UTD and accepts material properties as input 
data, 1t 1s regarded as a physical model. Moreover, among other calculated parameters, it 
implicitly extracts the PL and therefore 1s a path loss model. No information about the 


response delay 1s computed. NEC-BSC 1s described in the next chapter. 
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Wl. NEC-BSC 


A. SHORT DESCRIPTION AND BACKGROUND 


The Numerical Electromagnetic Code-Basic Scattering Code (NEC-BSC) is a 
high frequency EM analysis code, which employs the principles of UTD while interfacing 
concurrently with the MoM. It can simulate complicated scattering structures by uulizing 
relatively simple geometrical models like cylinders, plates, spheres etc. It was developed 
by ihe Ohio State UCniversiqtleciroscience Lab (OSU-ESI) in the carly 1980's umdes 
U.S. government contract and has been upgraded many times since then. 

It was originally wntten in FORTRAN 66 (version |) but all subsequent versions 
were created with the use of improved FORTRAN 77. Very recently, (September 1999) 
FORTRAN 90 conversion commenced [38]. It is considered to be a user-fnendly 
program since it requires only an elementary theoretical background in the fields of UTD, 
MoM, and EM diffraction. While its main features are presented later in this chapter, we 
should clarify at this point that this 1s not an antenna code in the sense that it does not 
culculate the current distribution or radiation resistance or any other pertinent information 
of antennas. Instead, given a specific current distribution, emphasis 1s given on the 


antenna pattern in various scattering environments [25]. 


b. OTHER WORKS BASED ON NEC-BSC 


So far. the code has been used with satisfactory results for the calculation of the 
total electric field and the near and far field antenna pattern that arises from the following 


antenna-scattering Structure interaction situations [39]. 


Zo 


e From the mast or other superstructures of a ship [40]. 

e From atank’s body 

e From the fuselage/wings of an aircraft 

e From the living quarters of the International Space Station 

NEC-BSC has also been used for comparison of measured and calculated patterns 

for different antenna placements and polarizations in these cases: 

e International Space Station [41] 

e Designing antennas for vehicular cellular applications (e.g. a patch- 
antenna studied tn [42]) 

e Development of the Advanced Enclosed Mast/Sensor system (AEM/S) 
project of U.S. Navy, first introduced on USS Arthur W. Radford in 
September 1997 [41]. 

e H-60 Helicopter antenna [43] 

e RCS measurements [44] 

A research in the literature has revealed that this code was not thought to be 
altracuve for predicting radio propagation in urban areas, with only one exception: a 
project in progress at OSU-ESL on cellular phone propagation around a building. Indeed, 
the code’s author asserts that this 1s the case as there 1s no other work on indoor radio 
propagation using NEC-BSC version 4 [45]. Under this perspective the approach was 
challenging. Version 4 accompanied by a Graphical User Interface (GUI) named NEC- 


BSC Workbench was available dunng this thesis work. It should be stressed, however, 
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that during this work throughout 2000 the code was upgraded at least three times 


illustrating in a vivid manner its continuing revision by OSU-ESL. 


C. NEC-BSC CAPABILITIES AND INHERENT ADVANTAGES 


NEC-BSC can be run either on an IBM-compatible PC using Windows as 
Operating system, or on a UNIX system. The GUI available for this code is the NEC-BSC 
Workbench and a Navy developed Silicon Graphics GUI for PC and UNIX, respectively. 
The exchange of data 1s performed in the form of ASCII files. The user inputs the data 
using either the Windows WordPad or Notepad, or the Workbench. During this process, 
which 1s based on a command word system, the user defines the geometry of structures 
and antennas and determines the various parameters needed to obtain a particular output 
antenna pattern. Initially, through the Workbench, which acts as a helpful editor, the user 
1s informed about whether or not a syntax or another similar error exists 1n the input file. 
Then, the program’s calculation engine is launched and generates warning messages, if 
any. and a progress bar. Finally, an output file will be created including the calculated 
results along with other files that the user might have asked for (e.g. an ORY file which 
permits the graphical representation of the EM rays used in the simulation). For instance, 
if the near-zone pattern of an antenna in a complex scattering environment was pursued. 
the resulting output will be the magnitude, phase, and dB value of the total electric and 
magnetic field, and the dB value of the Poynting vector at predefined observation points. 

The Workbench features a visualization too] which greatly facilitates the geometry 
building. More specifically, the user builds the complicated environment by using simple 


shapes (plates, cylinders etc) and checks the arrangement with the help of Workbench. 


\) 
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Once the simulation run 1s over. the same graphical interface can be used for visualizing 
the reflected/scattered/ditfracted paths. depicted as EM “rays”, and various plots (e.g. 
total electric tield versus distance). 

The most significant capabilities of NEC-BSC are [39]: 

e It can predict the Far/Near field pattern of an antenna in a complex 
environment 

e lt can provide the EMC or EM coupling between antennas in various 
Situations 

e [tcan determine potential radiation hazards 

e [lt can optumize the design and placement of antennas at the early stages of 
development 

e [t can furnish the user with Radar Cross Section (RCS) calculations 

e [t includes the option for single frequency or multiple frequencies 
simulations 

e UTD single, multiple, and triple plate interactions are possible. 

The above discussion should enable the reader to understand the majority of the 
inherent advantages of NEC-BSC. However, the rest of them become clear only with 
experience. A good example in this respect is given in [42], though, we should note, the 
authors do not allude to the version they implemented. In any case, it would be instructive 
to consider their deductions as regards NEC-BSC: 


e [tis fast and time-efficient 
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e Various vehicle elements can be easily added or removed allowing for 
ditferent EM considerations 


e The position of the sources can be arbitrarily changed 


D. LIMITATIONS OF THE CODE 


Many of the code’s Jimitations are directly caused by the very nature of UTD, 
Mamiels 2a |: 

e The scattering mechanisms need to be known and included in the model. 
Should some be excluded or accidentally left out, the accuracy of 
modeling will be gravely affected. 

e UTD may be employed only when the elements making up the model are 
large in terms of the wavelength. 

Furthermore, according to NEC-BSC user’s manual [39]: 

e Each plate should have edges at least a wavelength long. For curved 
surfaces the requirement is that their major and minor radii and length 
should be at least a wavelength in extent. 

e The source must be at least a wavelength from the surface and each 
antenna element should be at least a wavelength from all edges. Also, the 
receiving elements cannot be placed on a plate. 

e The dielectric layers must be thin and support only one mode of the fields. 

e At present, the solution does not contain surface waves, and hence the 


antenna cannot be placed close to the surface. 


LO 
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e The code 1s inadequate if a far zone backscatter or bistatic scattering result 
1s desired. 
e <A diffracted field from the plate-curved surface junction 1s not considered 
in this version 
e Accurate results can be expected for at most the first 30 dB of the pattern 
e Graphics generation through the Workbench 1s restricted to some 2D (two 
dimensional) plots only. For 3D (three dimensional), polar, and other more 
complicated plots the user should resort to exporting the numerical results 
to other applications (e.g. MathCAD. Matlab etc). 
e According to special instructions given in [38], there are certain limits as 
far as the total number allowed for each vanable is concerned. These 
limits are shown in Table 3.1. 
Numerous as these limitations might appear to be, they did not inhibit, at least 
severely, the indoor propagation simulations. All objects’ dimensions in the input file 
were much larger than the wavelength and this fact alone satisfied most requirements. 


Further discussion 1s made in Chapter V. 
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Table 3.1: Vanables limits 


















jk 


TAS PAGE INTENTIONAL YORE T Biba 


a2 


IV. NEC-BSC EXAMPLES 


This chapter presents three simple examples of indoor prediction using the NEC- 
BSC. The objective 1s to demonstrate some of the code’s capabilities in determining the 
neur-filed pattern of an antenna surrounded by obstructions much larger than the 
wavelength. First, the geometry 1s defined and a short description of shapes and materials 
comprising the model is made. Then, the observation points for each run are explained. 
Finally, the extracted results are shown in various plots and their interpretation 1s 
discussed. MathCAD was used for these and all subsequent calculations and plots shown 
in this thesis. In only one case. shown later in this chapter, the NEC-BSC Workbench was 
used and this was done for demonstration purposes. Al] input files used in the following 


simulations are given in Appendix A. 


A. A RECTANGULAR PEC ROOM WITH A CYLINDER 


In this first example, we calculate the electric field (E-field) for two different near- 
field patterns: one forming a circle around an obstacle and one running through the 
middle of the room. We consider a room with dimensions 6m x 14m x 6m (width x 
length x height) containing a cylinder of | m height and 1 m radius. The cylinder is 
located precisely at the center of the room. All room’s walls. the floor and ceiling, and the 
cylinder are perfect electnc conductor (PEC) materials. Midway between the cylinder’s 
central point and the 6-meter wall a vertically polamzed half-wave dipole is placed. The 
frequency of transmission 1s 2.45GHz and hence the wavelength is 12.2364cm. In order 


to contorm to all NEC-BSC requirements explained in Chapter 3, the dipole 1s positioned 
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at 3m (1.e.. 24.5 umes the wavelength) away from all obstacles. Furthermore, to account 
for correct attachment of the cvlinder to both the ceiling and floor we let the cylinder 


intersect the floor and the ceiling at depth 0.07 umes the wavelength [39]. [43]. 








Antenna 


Figure 4.1: Ground Plan of Room with Cylinder 

Ji Circular Near-Field Pattern 

The E-field in 360 points (each one corresponding to 1 degree step) around the 
cylinder is first calculated. The answer is given in dBV/m for all E-field components as a 
tuncuon of the observation points coordinates. A particularly helpful visualization tool of 
the Workbench allows the user to actually see the area covered by the rays launched by 
NEC-BSC (figure 4.2). This information is obtained by using the command LY and the 
program will generate an ASCII file by the extension ORY, which compnises the 
posiuons of the starting. hit, and end points of the rays along with signal strength values 


(39] It should be noted that by virtue of user defined filters 1t 1s possible to exclude all 


rays smaller than a certain value. Moreover. the number of interactions shown can be 


adjusted. 


Area covered 
by rays 





pee yy = 


Figure 4.2: Circular Pattern: Area Covered by Rays 


In figure 4.3, the plot of the results for the three E-field components as a function 
of angle g 1s shown. As expected due to symmetry, the values at angles 0 and 180 degrees 
are equal. However, it 1s observed that the signal strength 1s larger at points between the 
cylinder and the antenna (g=270 degrees) than behind the cylinder (g=90 degrees). This 
fact can be attributed to the cylinder acting as a strong reflector in the former and as an 
obstacle in the latter case. The Workbench creates a plot of all three components but not 
of the total E-field. Hence, a separate application should be used if other than the default 


plot 1s desirable. 
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Figure 4.3: E-tield Components versus Angle g Plot Created by Workbench 


pss Line near-field pattern 


Let us now find the variation of the E-field along a line passing through the 
centers of the antenna and cylinder. The observation points are sampled at a rate of half- 
wavelength (1.e.. 0.0612m): therefore, a line consisting of 229 points 1s formed. Again the 
treld in tront of the cylinder (where the antenna 1s located) is clearly stronger compared to 
the area behind the cylinder. In addition, the multipath phenomenon caused by the 
numerous reflections off the walls and the diffraction from the cylinder surface accounts 


tor the deep attenuations observed in the plot shown in figure 4.4. The discontinuity at the 
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7-meters distance point is due to the presence of the PEC cylinder. Ideally, no penetration 
of ravs inside the cylinder should occur. However, it can be noticed that there is erroneous 
prediction of E-field existence in !m-depth inside the cylinder’s interior, which manifests 
itself in the plot shown in figure 4.4. Every effort was taken to eradicate this anomaly. 
Despite the fact that the cylinder was slightly tilted and rotated in several consecutive test 
simulations. the problem persisted. No explanation could be found for this discrepancy 


upart from the suspicion that the code probably suffers from an undetected bug. 
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Figure 4.4: Total E-field versus Distance (vertical polarization) 
b. A TYPICAL ROOM 
In this example, a volumetric E-field pattern of an antenna located inside a typical 
room 1S Calculated. All room's dimensions and materials’ constitutive parameters were 


found in [47] and are shown in Tables 4.1 and 4.2. However, the observation points were 
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redetined to meet this example’s particular needs. The room along with the pattern used 
for this simulation 1s shown in figure 4.5. 

The antenna used is a vertically polarized dipole with piecewise sinusoidal current 
distribution and the transmission frequency is 9OOMHz, hence the wavelength is equal to 
0.3331m. Again no obstacle exists within the antenna’s near-field. This ttme we would 
like to know the excess gain/loss of the E-field in the entire room at a level 1.22m above 
the ground. For that purpose a volumetric pattern with the outer and inner loop consisting 
of 54 and 32 points, respectively, 1s defined. The sampling rate 1s half-wavelength which 
corresponds to approximately 0.1665m. Thus, the E-field is calculated at 1728 
observation points. The resulting “excess gain”, 1.e., the ratio of the total E-field in room 
over the total E-field in free space for every observation point, is shown in figure 4.6. It 
should be noted that for this and all subsequent excess gain/loss calculations a dynamic 
range of 30dB has been applied. In other words, all values more than 30dB lower than 
the maximum value have been set to 30 dB below the maximum. This practice was 


lollowed in order to comply with the NEC-BSC limitations [39]. 


All di | | 
inches 2 - 
room 


Table 4.1: Room’s Dimensions 


Sis 


Some different statistics were also considered in this room. First, the vanation of 
the normalized difference between the predicted and the free space E-field as a function 
of distance from the antenna was investigated. Accordingly, the whole volumetnic pattern 
was divided into several circular discs with radii ranging from 0.5 to 7 meters, thereby 
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creating 14 “circular sectors.” Hence, all observation points fall in one of these sectors 
forming various “groups.” Each of them corresponds to a certain distance from the 
untennu. An important fact is that not every group has the same number of points owing 
to the unequal surface of each circular sector. Nevertheless, very significant conclusions 
can be reached by this approach. In figure 4.6, the resulting histogram of the 
aforementioned normalized differences 1s shown. The X-axis is the radius of each circular 


disc trom the untenna in “steps” of 0.5 meters while the Y-axis is the number of “bins” 


tor this computauon. 


Material | Thickness Relative Dielectric Relative Magnetic 
(in inches) | permittivity loss permeability 


2"° layer: 
6°” bricks 


rd 


3” layer: 
G6 bricks 





Table 4.2: Room’s Properties 


ayy, 


Second, the probability density function (pdf) of the normalized difference 
between the E-field in the room and in free space for all distances 1s estimated and plotted 
ugainst a Gaussian pdf for direct comparison. The Gaussian distribution, shown in figure 
4.7. has zero mean and 0.6 times the variance of the normalized difference random 


variable. 
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Figure 4.5: A Typical Room 
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Figure 4.6: Excess Gain (dB) (oblique view) 
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Figure +.6: Excess Gain (dB) (top view) 
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Figure 4.8: Normalized Difference Random Variable Compared to Gaussian pdf 


Finally, the cumulative distribution function (cdf) for the normalized difference 
random variable for all ranges as a function of the loss in dB 1s determined. This is useful 
for communicauion engineering purposes and it is shown in figure 4.9. In the same figure 
ican be nouced how closely the Gaussian cdf matches the normalized difference random 
variable cdf. Further discussion of the importance of these statistics is made in Chapter 
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Figure 4.9: Normalized Difference Random Variable cdf versus Loss in dB 
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c. COMPARISON WITH THE SHOOTING-AND-BOUNCING RAY (SBR) 
METHOD 


In this example, the path loss along a line inside a rectangular room containing a 
closet 1s found. The dimensions and material properties of this room are described in [48] 
where the path loss was extracted using a modified shooting-and-bouncing ray/image ray- 
tracing approach (SBR). The same room was also modeled in [31] and [37] and the 
results were compared to those of [48]. The antenna used was a vertical Hertzian dipole 
and the operating frequency 1!GHz. A visual representation of the NEC-BSC model is 


shown in figure 4.10. 
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Figure +.10: Room with Closet 


As 1t 18 usually the case for most deterministic models, the task of choosing the 
suitable matenal parameters for accurately modeling composite or even uniform walls 1s 


challenging. More specifically, in NEC-BSC special consideration is required in selecting 
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the appropriate ZD commands as well as the correct relative permittivity (€,), loss 
tangent (tano ), and relative permeability (,) values for all walls and obstacles, the 


floor and the ceiling. This fact has been venfied many times during these comparison 
simulations. Every effort was taken to detect that combination of commands and material 
properties, which would lead to the closest possible match between NEC-BSC and SBR 
method’s predicted path loss. Clearly, since the NEC-BSC user’s manual does not but 
scarcely lend itself into indoor environments modeling, all the burden falls on the user 
who will have to have a good deal of pnor EM indoor propagation exposure. 
Furthermore. for the same reason, the correct interpretation of the results and/or the 
efficient intervention in the input file to make alterations, wherever applicable, are 
sufficient prerequisites to exclude inexperienced users form working on NEC-BSC. Prior 
(oO ulriving at the results shown tn this section, many discouraging predicted data were 
found und extensive investigation had to be performed in order to uncover the cause of 
the problems even for such a relatively simple room as this one. 

In the second example, the floor was modeled as “half-space” material [39]. In 
this case, however. this choice did not work. Moreover, the information provided in [48] 
strongly indicated that the use of “one-sided coated” material ZD-command for modeling 
the ceiling and floor was necessary. Again this approach was not justified by the results. 
Also, converting all plates into PEC material type yielded rather optimistic results for the 
first several meters and inaccurate results thereafter. Therefore, contrary to what was 
originally perceived as the mght path to follow, the employment of transparent slabs for 


all room's plates was chosen, resulting in satisfactory agreement with the SBR method. 
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Yet this approach assumes that no PEC backing behind the slabs is present, which 
could be viewed as a shortcoming intentionally imposed on the model. To lend credence 
to this unorthodox ussumpuion. two different input files were created: one with PEC slabs 
ul distance 0.00667 times the wavelength (0.299792 meters) behind the walls and another 
without these PEC slabs. The idea was to try an indirect method of accounting for the 
PEC backing by inserting 1n the original input file PEC slabs, allowing at the same time a 
smal] vet non-negligible “air-cushion” in between the wall slabs. The resulting path loss 
of the two approaches 1s plotted against the SBR results for direct comparison (figures 
4.11 and 4.12). Furthermore, an even more illuminating plot is given in figure 4.13 where 
the relative accuracy, meaning the difference (in linear values) between the path loss 
values predicted by BSC and SBR, of the PEC and no-PEC results are contrasted with 
euch other. The sharp 1:mprovement achieved by the “no-PEC solution” manifests itself by 
the relative proximity of the “no-PEC curve’ to the x-axis (distance) as opposed to the 


lurger variations over and below the x-axis of the “PEC solution” curve. 
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Figure 4.11: SBR and NEC-BSC (w/o PEC) Results 
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Figure 4.12: SBR and NEC-BSC (with PEC) Results 
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Figure 4.13: Contrasting the Relative Difference of the PEC with the no-PEC 


Solution 


Having found a particular NEC-BSC model. which more accurately described the 
path loss variation along a line inside the furnished room, the next course of action was to 
use that model to draw some additional conclusions with regard to both NEC-BSC 
validity and indoor propagation in the room under consideration. In doing so, the figure 
4+ 14 1s drawn where the path loss versus distance for the empty (denoted Plempty) and 
the furnished (denoted Pifurnished) room are shown. The difference between the two can 
be explained by the existence of a typical multipath-fading environment being further 
affected by the presence or absence of the closet for the furnished and empty room, 


respectively. Furthermore. the path loss in the furnished room 1s slightly less than the 
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path loss tn the empty room for all points before the closet’s location (at z = 11m). By 
contrust. the former 1s significantly larger than the latter right after the furniture. As in the 
first example, this is due to the fact that the closet is a strong reflector for the part of the 
room closer to the antenna whereas it 1s a blocking obstacle for the remainder of the 
room, 

Finally. in figure 4.11. the NEC-BSC predicted path loss 1s compured to the SBR 
method's results where a relatively good agreement is noticed. It should be noted that 
higher sampling rate was used in NEC-BSC compared to SBR depicted data. In 
particular, in NEC-BSC there was | observation point every lcm while | data point every 
20cm was used in drawing the SBR curve. In addition, the SBR data points were taken 
directly from the path loss graph of [48] rather than from a raw data file and, 
consequently, errors are certain to have been introduced in the approximations applied 
during that process. Therefore, some differences observed in figure 4.11 can be partly 
attributed to this fact. In any event. even in the comparisons presented in [31] and [37] 
there was not exact match between the curves of the other methods and the SBR solution. 
In conclusion, though not impeccable, the NEC-BSC results obtained under the afore- 
described assumption of no-PEC backing are reckoned to be acceptable adding some 


confidence to employing this code in other indoor environments. 
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Figure 4.14: NEC-BSC Path Loss Results for the Furnished and Empty Room 


¥. THE MODEL 


No simulation or test 1s complete without a description of the environment under 
examination. Equally essential is the inclusion of all indispensable simplifications that 
had to be made in order to conduct the evaluation of the too] in question. These topics are 


covered in this chapter. 


A. GEOMETRY DEFINITION 


The measurements reported in [11] took place on the second floor of Bullard Hall 
at the Naval Postgraduate School in Monterey, California. This building 1s rather old, 
built in 1958. The second floor of this building has enough complexity to create a harsh 
multipath environment for indoor radio propagation. Its dimensions are 39.8272m x 
39.8272m x 3.81m (length x width x height) and it contains a vast diversity of materials 
like concrete for exterior walls, glass for windows, plasterboard for most partitions. 
sheetrock with metal studs for some walls, acoustic tiles, wooden tables, concrete 
reinforced columns and more. [ts main function is providing space for laboratory 
experiments, therefore a large amount of testing equipment is spread around its main hall, 
which occupies a large part of the total available space. However, it also includes several 
offices and classrooms (with desks, chairs, wall-shelves etc), which are located at the 
perimeter leaving the laboratory at the center of the second floor. Hence, both lossy and 
good retlector objects are present throughout the second floor of Bullard Hall posing a 
profound problem for all modeling tools. A visual representation of the Bullard’s second 


tloor created by NEC-BSC Workbench 1s shown in figure 5.1. 


Antenna location 








Figure 5.1: Bullard’s Second Floor Created with the NEC-BSC 


Defining the geometry of the NEC-BSC model in this particular problem was not 
an easy task. Because of the age of the building CAD data was not available. Instead, the 
building s blueprints were used to extract al] necessary information about dimensions and 
materials inside Bullard Hall. That data was far from sufficient, and on-site inspection 
Was made in order to ensure the best possible accuracy. Yet it was not possible to model 
some significant construction details and thus a number of compromises were made, 
which are described later in this chapter. It should be emphasized that even if data in 
CAD format were readily available, 11 would not have been of much help. The current 
NEC-BSC version is not accompanied by CAD compatible graphical user interface 
(GUI), which would greatly facilitate the ASCII data input. It is assumed that a drawing 
tool developed particularly for NEC-BSC would minimize the time required for creating 


the mode] [46]. 
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Once all dimensions and coordinates were idenufied, the Bullard’s second floor 
was modeled by employing the single or multiple plates commands. That way, all exterior 
walls and windows, the floor and ceiling, all office and room partitions, all wooden 
tables, and most columns and sheetrock walls were included in the model. In wnting the 
Input file. special attention was paid to the area near the antenna whereas most details 
were omitted further away: they were assumed to have an insignificant effect on the radio 
propagation at close distances from the antenna. Furthermore, in developing the geometry 
of the simulated floor. one exceedingly favourable fact, which was taken advantage of, 
was the existing symmetry of many structures. For example, the dimensions of all 
office/clussroom partitions were identical. In such cases, the employment of the RT 
command expedited the entre work. In short, by virtue of this command, the user has to 
create only once the plate or structure which repeats itself several times. Once this 1s 
done, the object can be duplicated as many times as it is necessary and shifted to the 
desired location. However, since this capability was not documented in the NEC-BSC 
manual, several tests had to be performed to confirm its validity. These tests were 
successful and proved that this “trick” yields correct results. In a different arrangement, 
for instance, it can be used to create the steps of a staircase faster than using geometry 
commands repeatedly. 

The antenna used in the model was a Hertzian dipole with uniform current 
distnibution and operating frequency 2.45GHz. It was located at the center of the room 
Where the actual access point was during the measurements of [11]. Although the 


polarizauion of the antenna is not reported in [11] it is assumed to have been vertical since 
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itis well Known by many experments that vertical polarization has an advantage of at 
least 4 dB tor same-floor indoor radio propagation [50]. Also, the antenna’s height above 
the ground was assumed to be | meter. In fact. in measurements reported in [3]. [16], 
[29], and [31] the antenna’s height ranged from | up to 3 meters, thereby giving credence 
to that assumpuon. 

It has been argued countless times in the literature that the correct selection of the 
consutuuve parameters of the walls, floor, ceiling and obstacles plays the most significant 
role in all physical models. NEC-BSC is no exception to this well-established rule. 
Unfortunately, the exact values for every single construction matenal are not available. 
The immense diversity of both composite and non-composite materials found in all 
buildings makes it almost impossible to tabulate them all. Hence, it 1s most of the times 
the cuse to approximate these values using various methods. In this model, the standard 
values for the most common materials used (cement, plaster. glass etc) were taken either 
from [49] or [29] and they are shown in table 5.1. In all cases the relative permeability 


and the magnetic loss tangent were assumed to be equal to !.0 and 0, respectively. 


(Jt 
o> 


Dielectric loss tangent 


Material Relative permittivity (€, ) 





Caneiete 0.004363 

Reinforced concrete 5.84 

Plaster ae 
Aruficial composite wall 0.004363 


Table 5.1: Constitutive Parameters Used in the NEC-BSC Input Files 
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Moreover, 1n two particular cases a different concept was introduced: in front of 
the instruments room where a metallic mesh is located, and in one wall of the antenna’s 
room, which consists of approximately 80% plaster and 20% metallic beams. The goal 


was to calculate the €. 1n these cases under the assumption that they be replaced by 


“artificial walls” with the same EM reflectance and transmittance and proportions of 


constituent materials. This concept is discussed next. 


1. Metallic Mesh 


On-site inspection of the metallic mesh in front of the instruments room revealed 
that if A and a are considered being the total surface of the mesh and the metallic 


structure, respectively, then a@/A equals 0.17355. This conclusion justifies the assumption 


a 
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that rf the actual metallic mesh was to be replaced by an “equivalent” wall having the 
sume material proporuons then the latter would consist of 82.6447% air and 17.3553% 
metal. It can be further assumed then that only the “‘air-part” of this so-called “equivalent” 
wall contributes to transmission and only the “metallic-part” reflects back EM energy. 


Therefore. if P| and P. are the incident and the transmitted, respectively, power of an 


ir 


EM wave impinging normally on the “equivalent wall,” and if no absorption occurs, then: 





P =(1-—)*P *(1-(F_|°)] (5.1) 
A 
P =P, *[1-(Does I (5.2) 


where [and Tare the reflection coefficients of air and mesh, respectively. It 


air HES 


ts also true that: 


ia = | _ \V Cran - 3) 
mesh 
l + V Sask 


= (5.4) 


where €), is the relative permittivity of the “equivalent wall” replacing the 


mesh. 


From equations (5.1) to (5.4) it follows that: 
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E ) =2-096, which is the value that will be used in the input file for the NEC- 


BSC Bullard Hall model. 

DE Equivalent’ Composite Wall 

Following the same concept of “replacing” the actual wall by an “equivalent” wall 
with a supposedly equivalent effect on the transmitted EM energy, the sidewall next to 


the antenna was assumed to consist of 80% plaster and 20% metallic studs. Consequently, 


the following equations can be derived: 








B SUS02 >> SU (P alll (5.6) 

P, =P, *(1—(Pimp| I (5.7) 
| = E REST! 

_ Lev Ente i (5.8) 
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Eek 
lee = 1+ JEimp (5.9) 


and [- are the transmission coefficients of the “plaster-part” of 


plaster COMP 


where I 
the actual wall and the “equivalent” composite wall, respectively. Also, €,,,,,., 18 the 
relative permittivity of plaster equal to 2.3, and €,,,,, 1s the relative permittivity of the 


“equivalent” wall with the same EM properties as the real one. From equations (5.6) to 
(5.9) it follows that: 
Cm acs 47 (5.10) 


qedetp 


and that was again the value used in the simulations pertaining to Bullard Hall. 
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B. MEASUREMENTS 


In 1999, laboratory testing was conducted on the 2" floor of Bullard Hall and the 
results of the signal strength and throughput measurements were reported in [11]. The 
goal was to evaluate four commercially available wireless networking components in a 
typical multipath environment. Therefore four different sets of data were obtained which 
give a good idea of the actual coverage area. 

One access point located in the thesis students’ room and two wireless client 
computers were used for the tests. The transmitted power ranged from 300 to 365mW 
depending on the type of component. As far as field strength results are concerned, they 
were measured either in percentage or as signal to noise ratio. The maximum coverage 
area was determined by detecting the locations at which the client was no longer able to 
maintain COmmunications with the access point. This fact, of course, limited the extent of 
measurements as a large portion of the second floor was ruled out and no data were 
obtained. Furthermore, the signal strength was not measured in the neighbouring to the 
antennas location rooms. The best performance was delivered by the Lucent 


Technologies WaveLAN IEEE components and the pertinent coverage is shown in figure 
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Figure 5.2: Lucent Technologies WaveLAN Coverage (from Ref.[11]) 


C. SIMPLIFICATIONS 


Many practical and code-related difficulties were encountered while developing 
the NEC-BSC model. Most of them were trivial and a way around was often found. 
However. some problems were persistent and almost impossible to overcome and a kind 
of compromise was made. In addition, in other cases for the sake of accelerating the 
simulations some simplifications were included in the model, which, strictly speaking, are 
likely to have affected the accuracy of the predicted results. In what follows, all the 


above actions are identified and descnbed concisely. 
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l. ZD commands 


Whenever an object inserted in the model has different properties (either 
consululrve parameters or thickness) than the one before, 11 must be preceded by an 
appropriate ZD command. Since the current version's configuration allows a maximum 
number of 25 ZD commands in each input file, the order of alternate geometry and 
Structure definiuon commands must be meticulously arranged so as to avoid error or 
warming messages from the code. Furthermore, being restricted for the above reasons, the 
user Will have to limit the overall number of different materials used in the model. For 
example, 1f within a particular room there are more than 25 different types of materials, 
the user will have to devise an “intermediate” set of properties for some of them to 
account for this limitauon of the code. 

eae Plates 

The model developed in this thesis comprised numerous plates. Various structures 
and furnishings encountered in Bullard Hall were excluded from the model, nonetheless, 
and one reason was the maximum allowed number of 200 plates for each input file. 
Another reason was the difficulty to attach together certain plates. It was found that the 
suggested plate attachment process described 1n [39] could not prevent the appearance of 
problems in some cases. Although all criteria were met in various plate attachments, the 
code sul] returned meaningless results sometimes, but not always accompanied by 
warming messages. This problem could be avoided by either shifting the plate in question 
Slighthh away from the originally desired position or by completely removing it from the 


input tile. Despite the above-described difticulues. the second floor NEC-BSC model did 
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not deviate from the standards of other applications (like the one in [29]) as far as 
seometrical representation 1s concerned. 


by Presence of people 


It has been shown that the density of people’s presence inside a building has a 
profound effect on the received signal statistics [51]. Even when the transmitter and 
receiver are fixed the human motion can result in 2-10 dB attenuation of the average 
signal strength [52]. This fact could not be taken into account. 


4. Radiation scattered by neighbouring buildings 


EM radiation can leave a building through a window and come back after being 
scattered by neighbouring buildings [33]. This fact can be crucial for larger distances as 
these “out-and-in” ray paths can carry substantially more power compared to rays 
propagating through many walls within the building. The neighbouring buildings were 
not included in the model, but their actual distance from Bullard Hall is fairly large 
(ranging from 30 to 60 meters) and therefore it is thought that their effect 1s negligible. 


>: Clear space boundaries 


In [32] it is argued that the boundaries of clear space, through which the signal 
propagates. plavs un important role and it should be included in any model. Whereas the 
lower bound was implicitly taken into consideration by including some furnishings, the 
Same was not true for the upper bound which is formed by pipes, support beams, 
ventilation ducts, and lighting fixtures that protrude from the ceiling of Bullard Hall's 


second floor. 
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VI. SIMULATION ANALYSIS 


A. PLATFORM CONFIGURATION AND TIME PERFORMANCE 


All simulations were executed on an IBM-compatible PC using Windows as its 
operating system. The CPU was an Intel Pentium II at 600 MHz and the RAM memory 
was 768MB. The NEC-BSC Workbench and the main executable NEC-BSC FORTRAN 
program versions used were the 4.1.37 and 4.2-06, respectively. These versions were not 
the most recent ones but, rather, the most stable. In general, the Workbench performed 
well even when large files (several Mbytes) were analyzed. More specifically, there were 
only two cases where the Workbench either crashed or displayed erroneous figures. First, 
when very large rays files (usually a couple Mbytes) were called the Workbench 
operation stalled and reboot was necessary to resume calculations. Second, sometimes 
when a relatively big object was placed next to a small one the Workbench responded by 
showing the former in a clearly shifted position away from the correct one. This problem 
appeared only in the oblique (default) view option. 

As far as time performance is concerned, it was found that it was proportional to 
the number of observation points, the complexity of the environment and the relative 
location of the antenna. In particular, the number of observation points varied from 696 
points, for the comparison file, to 26,244 points for the statistics files. In the latter files 
the observation points were taken at a sampling rate twice the wavelength. For all 
practical purposes, that was the best attainable sampling rate, as anything better would 


require unreasonable CPU time resources. The total number of plates in the input files 
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was ISI. while the number of materials used was 19. Furthermore, the antenna was 
placed at five different locations including the room from which all actual measurements 
were carned out. It was observed that the time required for the calculation of the E-field 
in the 26,.244-points file ranged from 1,536 to 2,293 CPU minutes, whereas only 145 
minutes were needed in the 696-points file. The conclusion arising from these findings 
plainly proves that t1me performance explicitly depends upon the total number of 


observauion points. 


b. COMPARISON WITH MEASUREMENTS 


Since the results with regard to signal strength in {1 1] are presented in the form of 
range of dBm values for different coverage areas throughout the Bullard’s second floor 
rather than for specific locations, direct comparison between measured and predicted 
values 1s not possible. Instead, a qualitative comparison can be accomplished by visually 
comparing the coverage obtained by the in-site testing and the corresponding predictions 
of NEC-BSC. Indeed. by contrasting figures 5.2 and 6.1 one can infer that there 1s a 
relatively good agreement at least in terms of overall coverage. Figure 6.1 shows the 
predicted excess gain/loss for the same floor. It 1s clear that, while the signal is strong 
inside the antenna room. it rapidly attenuates along the Y-axis where penetration through 
many walls occurs. On the contrary, across the antenna room (along the X-axis), where 
only a couple of walls exist, the attenuation is far less and the signal reaches a greater 


distance before being reduced to low levels. 
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Figure 6.1: Predicted Excess Gain/Loss in Bullard 2°* Floor (open doors) 
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A second comparison approach was also conducted. The figure 5.2 was enlarged 
and printed out and then the whole floor map was divided into small squares. The side of 
each of these squares corresponded to 1.32 meters distance. The areas not marked in 
figure 5.2 were excluded. The antenna room was also ignored because there was not any 
range of values given for this area. Then, working on the remainder of the floor, a data 
matri\ was extracted whose elements were a function of location and respective signal 
strength value in dBm. Next, an input file with an equivalent layout of observation points 
Was generated and run and the resulting E-field values were converted to receiving power 


in dBm. Finally, the two matrices were plotted together in figure 6.3. 
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Figure 6.3: Comparison of Measured and Predicted Data Matrices 
While comparing the two data matrices, the numerous inherent errors of the afore- 
discussed procedure have to be taken into account. Large portions of the floor were left 
out and the measured data matrix was by no means accurate. Rather, “gracious” judgment 


was often employed in determining the dBm value of certain points. On the other hand, 
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the observation points of the NEC-BSC input file had be specified so as to approximate 
the arrangement shown in figure 5.2. That, in turn, 1s believed to have caused some 
additional errors. Furthermore, it is well known that due to deep fading in the indoor 
propagation, the E-field can vary considerably at distances only a fraction of the 
wavelength in any direction. Thus, if some of the observation points used in the NEC- 
BSC input file fell into these “deep-fading” regions, the comparison would have been of 
littke merit. Whether this has indeed happened in this particular case cannot be positively 
verified Owing to lack of position-specific measured data. Nevertheless, it can still be 
seen in figure 6.3 that along the X-axis the coverage area is more or less the same. The 
difference in the Y-axis could be probably explained by inappropnate modeling of the 


walls close to the antenna s§ site. 


Os SIMULATIONS 
I. Testing Different Model Approaches 


The orginal input file was slightly modified in two different ways and, 
consequently. two new files were obtained. The objective was to test the effect of these 
modifications on the predicted gain/loss map. In particular, in the first test, a couple of 
walls in the proximity of the antenna were modeled as “double-sided coated” material 
plates as opposed to transparent slabs [39]. The thickness and constitutive parameters 
information was retained and only the ZD commands were altered. This was done in 
order to find the most accurate method of modeling the concrete walls containing metal 
studs in the antenna’s close neighborhood. This concept failed in bearing meaningful 


results as it is clearly shown in figure 6.4. No propagation appears to be possible beyond 
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these walls. Therefore. this result. although expected since there 1s a PEC in the middle 
coated on both sides of the wall. does not agree with the measured data. 

In the second test. the doors of all rooms in the antenna’s side of the floor were 
assumed to be closed and the input file was modified accordingly. The resulting gain/loss 
mup 1s shown in figure 6.2. No discernible difference with the orginal opened-doors file 
can be seen, which probably implies that this scenario was not accurately modeled, or that 


the doors have only a modest effect on the signal propagation. 
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Figure 6.4: Predicted Excess Gain/Loss Using Double-Sided Coated Walls 


2 Comparing with Free Space 

Apart from the results obtained from the original antenna position input file and 
shown in figure 6.1, four additional output data files were also computed, each one with 
the antenna at a different location inside the Bullard 2°" floor. While that was done for the 


extraction of some interesting statistics described in the following section, some helpful 
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excess gauin/loss maps of the entire floor were also obtained and one of them is shown in 
figure 6.5. In this figure, resulting from the input file with the antenna in the middle of the 
floor, it can be seen that the predicted coverage area is significantly larger compared to 
that of figure 6.1. This fact could most certainly be expected for it has been proved 
through various measurements reported in the literature that, in the majority of same-floor 
propagation cases, the best coverage 1s achieved by placing the antenna in the middle of 


the floor. 
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Figure 6.5: Predicted Excess Gain/Loss (antenna in the middle of the floor) 


D. STATISTICS 
l. Flow of Work 
Following the approach already discussed in Chapter IV, the four additional data 


files for the four different antenna placements were analyzed in the manner described in 


this section. Once more, the objective was to extract the statistics of the excess gain/loss 


ee 


of E-neld relative to free space in the entire floor for various antenna locations including 
the iniual one. All antenna locations are shown in figure 6.6, which was created by NEC- 
BSC Workbench. MathCAD was employed in all cases and the flow of work was the 
following: 

e Data input: The program read the two NEC-BSC data files, 1.e., the one 
resulting from the furnished and free space input file, respectively. Two 
\wo-column matrices were generated from each data file; the first column 
contained the distance of the observation point from the antenna while the 
second one contained the value of the E-field for that specific point. 

e Data grouping: The entire observation-points layout was grouped into 
several circular sectors in steps of 50cm length. Thus, every group 
comprised those observation points, which fell in its respective sector. 
Each resulting group had different number of elements, which formed a 
set of data for every distance increment from the antenna’s position. 

¢ Room gain random variable: A random vanable was introduced in the 


following fashion. In every so-called “‘distance-group” the ratio of the E- 


field in the room (denoted E 





) and the E-field in free space (denoted 


FOC 
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iE) was calculated by 
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: h 
where m and n denote the mm” distance group and n” particular value for 


this group, respectively, and KR, corresponds to the value of the pertinent 


Al 
“room gain” over free space. Consequently, every distance group created 
one realization (out of m in total) of the so-called E-field room gain 


random variable. 
Probability density function (pdf): Let us assume R, to be the random 


variable describing the variation of the normalized difference for the m” 


realization. Then the pdf for all R, can be computed and the resulting 


histogram for all realizations can be drawn. In figure 6.7, the histogram 
obtained from the output file with the antenna at the center of the room 1s 


shown and denoted by L,, . 


Average pdf: Let us assume r to be the random variable descnbing the 
variation of the function resulting from averaging the pdfs for all distance 
groups. The ensemble pdf average of all distance groups was computed 
and its Statistics were investigated. A Gaussian pdf having the same mean 
and variance as the “average-pdf’ was generated. Both pdfs were plotted 
versus excess gain in dB on the same graph and they are shown in figure 
6.8 for the antenna located at the center of the floor. An excellent 
agreement was observed and it was concluded that the excess room gain 
random variable r in dB has the lognormal distnbution. This deduction 1s 


further discussed and explained later in this chapter. 
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Figure 6.6: Antenna Locations Used in the Statistics Analysis Input Files 
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Figure 6.7: Histogram of All Distance-Groups pdfs (Ant. at Location 2) 
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Figure 6.8: Average-pdf of the Room Gain Plotted Against Gaussian pdf 


Cumulative distnbution function (cdf): The cdf of the random variable r 
was computed and plotted versus the room gain in dB along with a 
Gaussian cdf again with the same statistical characternstics (mean, 
variance), and it is shown in figure 6.9 for the antenna at the center of the 
floor. Once more, the agreement was remarkable yielding another 
manifestation of the fact that /ogr has a normal (Gaussian) distribution. 
The merits of such graphs as the one shown in figure 6.9 are quite 
noticeable for the communications engineer for it depicts practical 


information used for the determination of the link budget expected in the 


12 


enure room. In particular, in the example shown 1n this figure, the points 
at Which the two marked dashed lines of 0.9 and 0.5 intersecting the cdf 
curve indicate that the estimated probability of encountering less than 


11dB and less than 4dB excess gain in the Bullard’s second floor 1s 90% 
and 50%, respectively. 
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Figure 6.9: Cdf of Random Vanable r (antenna at location 2) 


a Results 


Revertuung to the staustics of the room excess gain random vaniable r, it should be 
noted that the lognormal distnmbution has often been found to describe the amplitude 


variauon of the indoor radio propagauion channel. as reported in [14]. However, in our 
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analysis, yet another tesumony to that fact has been given through a different perspective. 
Contrary to the approach discussed in [14] and elsewhere in literature where the statistics 
applies to different distances, 1n our case they apply to the same distance as defined by the 
so-called “distance-groups.”” More specifically, this concept is employed for arbitrary 
locations of the transmitter and the “receiver” (i.e., where the observation points are 
placed), as long as the distance between them is the same and the simulation 1s conducted 
on the same floor. A thorough, yet not guaranteed to include everything, research suggests 
that this ““same-distance” approach has not been reported before. In any event, the large 
number of reflections and the spatial inhomogeneity of the mobile channel seem to be 
sufficient in explaining this statistical behaviour. More specifically, if each reflection is 
assumed to contribute by a factor (reflection coefficient) to the total magnitude (a product 
of any factors) then if the logarithm of the product is taken it will yield a sum of the 
logarithms of the particular factors. Therefore, by applying the central limit theorem, it 
can be inferred that for a large number of reflections the resulting signal magnitude will 
resemble a random variable with normal distnbuuon, which 1s precisely what it is noted 
in figures 6.8 and 6.9. As far as the inhomogeneity involvement is concerned in this 
process, it has been proved to contribute to the transition from Rayleigh and Rician 
distribution in the outdoor radio channel to lognormal distribution in the indoor radio 


channel [14]. The values of the “room gain” parameters used for all antenna locations 


comparison can be found in table 6.1. 
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Table 6.1: Room Gain Parameters 


The resulting average-pdf curves for all, except the initial, antenna placements are 

shown in figure 6.10. They were plotted against the corresponding Gaussian distnbution, 
e., having the same mean and vanance. Despite the better accuracy of the results 
obtained when the antenna ts located at the center of the floor, the curves gained from the 
other locations can stil] be considered adequate for venfying the Gaussian charactenstics 
ol the random variable r. Likewise, the same principle holds for the cdfs of the rest of the 


antenna placements shown in figure 6.11. 
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Figure 6.10: Average-pdf Plotted Against Gaussian Distribution 
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Figure 6.11: Room-Gain cdf Plotted Against Gaussian Distribution 


a: Other problems 


For the benetit of future work on indoor propagation using NEC-BSC, several 


problems not directly related to the scope of this thesis will be reported in this section. 


They were discovered either during the early familiarization stages or during the actual 


development of the input files. Although many more than those presented here were 


found, it was reckoned more useful to include only those that have not yet been resolved. 
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a) Command ZD 11 


Even though the command ZD 11, which 1s used to define absorbing 
material properties, 1s descnbed in [39] as being included in the code, it does not 
cooperate with the NEC-BSC Workbench. Specifically, the window of the ZD command 
does not list this option and, furthermore, if the user inserts it manually and then double- 
clicks on it the application crashes. No other cause than a potential bug could be found. It 
was hoped that this command option could be used to model several absorbing materials 


euler Fall, 


b) Commands not supported by Workbench 


It was noticed that the command wizard window of the Workbench did not 
support the commands BF, FT, LY, Pl. TA, TD, TW, TY, VM. Moreover, the 
Workbench supports the command ZC, but its function remains unknown since It 1S not 


included in the Help file. 


C) Polar plot command 


A minor yet worth reporting problem was detected regarding the command 
PP that is used to generate either a polar or a rectangular plot. In particular, the problem 
les on the logical variable LPPREC. If the user chooses the polar plot option in the PP 
command window, the Workbench ought to respond by specifying the value “False” to 
the LPPREC variable according to the Help file. However, this does not occur and the 
Workbench gives the value “True” automatically, which is the opposite of the desired 


value. In addition, even if the user persists and corrects manually the LPPREC value back 


aS 


to “False.” the Workbench will again change it 


the PP command. 





VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis dealt with the indoor radio propagation channel and treated the 
problems caused by the complex characteristics of the latter. Several modeling methods 
used in this field were discussed and their pros and cons were presented. The major 
objective, however, was the evaluation of the Numerical Electromagnetic Code- Basic 
Scattering Code applicability in confined environments. In that respect, the NEC-BSC 
results were compared to those obtained by other applications as well as actual 
measurements conducted in the Bullard Hall of NPS in Monterey, California. In general 
terms the evaluation, although not complete, provided invaluable information regarding 
the code in question and strong indications of the validity of the code in modeling the 
indoor radio channel. In what follows, the most significant conclusions are summarized 


und some recommendations for further study are given. 


A. CONCLUSIONS AND UNRESOLVED ISSUES 
e Despite the lack of previous work in indoor radio propagation channel 
modeling based specifically on NEC-BSC, it has been proved that this 
code ugrees to u satisfactory degree with the results reported in [48] for 
the same arrangement examination. Furthermore, a qualitative comparison 
between NEC-BSC predicted E-field values and actual measured signal 
Strength for the same building yielded promising results, thus warranting 
further investigation in other buildings and different indoor radio 


propagation experiments. 


oul 


In simple arrangements, easily meeting the code requirements discussed 
in Chapter III. NEC-BSC simulation completion was fast. Nevertheless, 
that was not the case in the Bullard Hall model where the run times were 
excessively high and proportional to the complexity of the environment, 
the number of the observation points and the relative placement of the 
antenna. Because of that, the sampling rate of the observation points had 
to be reduced below the usually applied “half-wavelength” rule of thumb 
and this 1s believed to have affected the accuracy of the results. 

The NEC-BSC user's manual [39] seems not to have included all plate 
attachment cases for there were numerous examples encountered during 
the development of the Bullard Hall model where the suggested process 
of intersecting/attaching two plates was followed yet the plate attachment 
clearly failed prompting the code to launch either error or warning 
messages. Sometimes the relevant code messages were made to vanish by 
shifting the plate slightly away from the desired position, usually a 
fraction of the wavelength. Whenever a way around this problem could 
not be found. the plate under consideration had to be left out of the model 
degrading that way the accuracy of the actual Bullard Hall representation. 

The E-field amplitude variation in confined environments has been 
reported to obey the Gaussian distribution statistics [14] when studied in 
terms of varying distance. The approach described in the previous chapter 


revealed that this 1s also true for the so-called “room excess gain” random 
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variable when studied in terms of same “distance groups.’ This fact, apart 
trom adding credence to the validity of NEC-BSC predicted values, 
implies that the lognormal distbution may well be applied in both 
varving and same distance statistics considerations. However, final 
conclusions on this matter can be reached only after further study 


(different rooms, materials etc). 


B. RECOMMENDATIONS 


Since the cuient NEC-55C€ craphical, ser interiace (GUI) ere the 
Workbench, does not support any CAD package, it 1s suggested that a 
future version of NEC-BSC incorporate AutoCAD interoperability, 
thereby enhancing the code’s functionality and _ user-fnendliness. 
Moreover. it 1S suggested that a database be built containing the 
consutuulve parameters of various frequently encountered walls and 
materials of the indoor radio channel. This database, if included in a 
future version of NEC-BSC, could speed-up the geometry development. 
The current NEC-BSC version does not allow insertion of all available 
commands through the Workbench. It 1s recommended that this capability 
be included in one of the future versions of NEC-BSC. 

A parallel execution option [37] could also be made available in an 
improved version of NEC-BSC for it would significantly shorten the run 
times. espectally in cases like the Bullard Hall model where the 
complexity of the environment ts critically high approaching the 
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limitations of the code illustrated in chapter III. Nevertheless, it should be 
said that this can be implicitly, but not so efficiently, performed by 
dividing up the computation “burden” of the observation points into 
different PCs (e.g. points 0-100 assigned on PC-1. points 100-200 
assigned on PC-2. etc). 

Whenever the near-field pattern 1s computed by NEC-BSC, the results 
include both the Electric and Magnetic field components along with the 
Poynting vector components. In case that only the signal strength is 
desired the computation of the H-field and the Poynting vector could be 
omitted in order to reduce the simulation time. However, this is not 
possible in the current version of the code and the user cannot choose the 
particular data to be calculated. 

As innovative technologies like Bluetooth [5] appear in the forefront of 
indoor applicauons and call for transition from nanocells to picocells, 
smaller room statistics will become necessary. NEC-BSC could be used in 
that respect. Furthermore. NEC-BSC could also be employed in surveying 
the room-to-room propagation and hallway-to-room statstics, which have 
not been taken into account in this thesis. 

In a future study, NEC-BSC could be utilized for extracting various co- 
polarization statistics in the same floor or room; these could include but 
not be limited to horizontal-to-horizontal, vertical-to-vertical and circular- 


lo-circular wansmitter-to-receiver polarization. Then their pertinent results 
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could be contrasted against each other and confirmed through 
measurements. In the same sense and for the same purpose the resulting 
cross-polarization (e.g. like vertical-to horizontal, horizontal-to-circular, 


etc) statistics could be investigated with the help of NEC-BSC. 
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APPENDIX A. NEC-BSC INPUT FILES 


CE: This is the input file used in the cylinder 
CE : room example 

CE:A rectangular room with PEC walls 
CE: containing a PEC cylinder in the 
CE: middle of the room 

eke 

2.45000 

UN: 

| 

LS 

0 

me Perec Electric Conductor 

0 

Gr Floor 

0 

PM: 

8 

0.000000. 0.000000, 0.000000 
O0.000000, 0.000000, 6.00000 
6.00000, 0.000000, 6.00000 
6.00000, 0.000000, 0.000000 
6.00000, 14.0000, 6.00000 

6.00000, 14.0000, 0.000000 
0.000000, 14.0000, 6.00000 
0.000000, 14.0000. 0.000000 

() 


Y) £& © 2 Oo) 


CoN Sse ey rare 


On @ oO -—) Chi 
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© wwn~) Uw Cc —1y ~) 


0 

CG, PEC. winder 

3.00000. 7.00000, 3.00000 

0.000000, 0.000000. 90.0000, 0.000000 
1.00000, 1.00000 

-3.00856, 90.0000, 3.00856, 90.0000 
Scr 

3.00000, 3.00000, 3.00000 

0.000000, 0.000000, 90.0000, 0.000000 
-1. 0.500000. 0.000000 

| OOOOO. 0.000000 

CE: For honzontal polarization 

CE: simply comment selection out 

CE: the above antenna settings 

CE: and uncomment selection out 

CE: the settings given below 

GM SG: 

CM: 3.00000, 3.00000. 3.00000 

CN\1:. 90.0000. 0.000000. 90.0000, 0.000000 
CMI: -1. 0.500000. 0.000000 

CE: 1.00000. 0.000000 

PN: Circular pattern around cylinder 
3.00000, 7.00000, 3.00000 

0.000000, 0.000000, 90.0000, 0.000000 
F 

2.00000, 90.0000, 0.000000 

0.O00000. 0.000000. 1 .OO000 

360 

CE: For line pattern simply comment selection out 
CE: the above PN command and 

CE: uncomment selection out the 

CE: following PN command 

CM: PN: Line pattern running through Cylinder’s Center of Mass 
CM: 0.000000, 0.000000, 0.000000 
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CM. 0.000000. 0.000000. 90.0000, 0.000000 
CMaai 

CM: 3.00000. 0.000000, 3.00000 

CM: 0.000000. 0.06 1 2000, 0.000000 

GE, 229 


ee 

i 

ae 

7 

OR 

IeIN: 

CE: End of cylinder room input file 
Glee 

CE: 

GE: 

Giz: 

CE: 

CE: This is the input file used in the 
CE: typical room example 

Gis: 

FR 

0.9 

UN 

3 

Ril: 

- 120.000, 0.000000. 0.000000 
0.000000. 0.000000. 90.0000. 0.000000 
VN: 

0.000000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
T 

6.55715, 6.55715, 48.0000 

6.55715, 6.55715, 0.000000 

1, 54,32 

Kh 

120.000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
ae: 

I 

Lay 

Er 

CMEaYy 


Sy 


CES. 

CMEPR 

CAS 

CE 5054557501 s 
US 

0 

5G 

O02 7 2.a.S. 
0..0..90..0. 

= OU 

| ..0. 

CM: RG 

CM: 0.,0.,0. 
CM: 0.,0.,90.,0. 
CI =? O2a20.0 


ae, 32000) li) 
PM CEILING 
16 
O05 -O.Ugrle.0 
20 5020 ie © 
See. sale 
EGO a OR Ole We 
120.0, 144.0, 112.0 
0.0. 144.0, 112.0 
44.0, 144.0, 84.0 
8.0, 144.0, 84.0 
202 020, a0 
44.0. 0.0, 14.0 
O05 2007-300) 
120,02 G02 050 
120.0. 144.0. 0.0 
44.0. 144.0, 0.0 
8.0. 144.0, 0.0 
0.0, 144.0. 0.0 


90 


= 


6 
2 
+ 
3 
2 
0 
0 
Ors BRICK 
| 
5 
0.53,3..0.048,1.,0 
Sol Ee SOR ee 
G3 ..0 048.1 ..0 
PM WALLS 
16 
OO MNOT OR) 223 8, 
e001 270 


en O01 12.0 
01050-0511 1220 
120.0. 144.0, 112.0 
0.0. 144.0, 112.0 
44.0, 144.0, 84.0 
8.0, 144.0, 84.0 
12.0, 0.0, 14.0 
44.0. 0.0, 14.0 
COLO On 0.0 
AOR OME ORO MEE ORO) 
120-0 es Or 00) 
44.0, 144.0, 0.0 
8.0, 144.0. 0.0 
0.0, 144,0, 0.0 


calle 


arp ees eer 


—=— HN CO = — JDO me RK ON Ce mt 
rel al A Eo A 


QO — 
CON — 


0) 

BM HALLS 

Ie 

120.0. U0) 0 

240.0, 0.0, 112.0 

240.0, 216.0, 112.0 

UO oO 1120 

=O 90102 be 
0.0. 0.0, 112.0 
120 0 Ow 
24050, (O00 800 

240.0. 216.0, 0.0 

OC, 2 ou 


-120.0, 0.0, 0.0 
ChCP CRO POR 
0 
| 
2 
§ 
0 
2 
3 
9 


7° GEASS 

| 

| 

Wale). 290013210. 


PM WINDOW 
16 
CHO ON OE Ja 


[70 020, 112:0 
44.0, 0.0, 112.0 
120.0. 0.0, 112.0 
120.0. 144.0, 112.0 
0.0. 144.0. 112.0 
44.0. 144.0. $4.0 
8.0. 144.0. 84.0 
P2000 1420 
44.0, 0.0, 14.0 
CO OU 020) 
ZN, ROH ONG 
120.0, 144.0, 0.0 
44.0. 144.0. 0.0 
§.0. 144.0. 0.0 
0.0, 144.0, 0.0 


ar? al i eee. 


Zhe 


Ap 

EN 

CE: End of typical room example 

CE 

Cr 

CE 

CE 

Che. 

Ck: This is the input file used in the 
CE: comparison with the SBR method 
Ce 

LN 

| 

es: 

0 

aK: 

| OOOOO 

TACs 

| 

0.130000, 4.00000, 0.0900000, 1.00000, 0.000000 
GEatioor 

0 

ZD: walls 

| 

2 

0.149800. 6.00000, 0.0300000. 1.00000, 0.000000 
0.100100, 6.00000, 0.0300000, 1.00000, 0.000000 
PM: 

8 

0.000000, 0.000000, 0.000000 

10.0000. 0.000000. 0.000000 

10.0000, 3.00000. 0.000000 

0.Q00000. 3.00000. 0.000000 

0.000000, 3.00000. 21.0000 

0.000000, 0.000000, 21.0000 

10.0000, 0.000000, 21.0000 

10.0000, 3.00000, 21.0000 

0 
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‘~~ lo_— 


POM he NODA OO ~)1 OW ~) OO WO L 


ZD: ceiling 

| 

| 

0.130000, 4.00000, 0.0900000, 1.00000, 0.000000 
OE JG; 

CM: 4, 0 

CM: 0.000000, 0.000000. 0.000000 

CM: 0.000000, 0.000000, 21.0000 

CM: 10.0000, 0.000000, 21.0000 

CE: 10.0000, 0.000000, 0.000000 

PG: 

4,0 

0.000000, 3.00000, 0.000000 

0.000000, 3.00000, 21.0000 

10.0000, 3.00000, 21.0000 

10.0000, 3.00000, 0.000000 

CE: For the “empty room" solution 

CE: simply comment selection out everything 
CE: until (but not including) the SG command 
ZD: closet 

QO 

RR: 


oS. 


1.50000, 0.000000, | 1.0000 
0.000000, 0.000000, 90.0000, 0.000000 
PM 

S 

0.000000, 0.000000. 0.000000 
7.00000, 0.000000, 0.000000 
7.00000, 2.00000. 0.000000 
0.000000, 2.00000. 0.000000 
0.000000. 2.00000. 0.600000 
0.000000. 0.000000. 0.600000 
7.00000, 0.000000. 0.600000 
7.00000. 2.00000. 0.600000 

¢ 


S&H CO Hh UW lv ~ 


43 OO WN 


om 


— ty Orv ~1 OO W 


~~] —/ 


Co LAN a —" GN SO). Ch Coe Ic 
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0 

RR: 

- 1.50000, 0.000000, - 11.0000 

0.000000, 0.000000. 90.0000, 0.000000 
CE: For the PEC backing solution 

CE: simply uncomment selection out everything 
CE: until (but not including) the SG command 
CM: RR: 

CM: 0.000000, 0.0020000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 
CM: PG: ceiling backing 

CM: 4, 0 

CM: 10.0000. 3.00000, 0.000000 

CM: 0.000000, 3.00000, 0.000000 

CM: 0.000000, 3.00000. 21.0000 

CM: 10.0000. 3.00000, 21.0000 

CIMA IRRe 

CM: 0.000000, -0.0040000, 0.000000 

CM: 0.000000. 0.000000, 90.0000, 0.000000 
CM: PG: floor backing 

CM: 4, 0 

CM: 10.0000. 0.000000, 0.000000 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000. 0.000000, 21.0000 

CM: 10.0000, 0.000000, 21.0000 

CM: RR: 

CM: 0.000000, 0.002000, 0.000000 

CE: 0.000000, 0.000000, 90.0000, 0.000000 
Gr. 

SG: 

2.00000, 2.50000, 3.00000 

90.0000, 90.0000. 90.0000, 0.000000 

-1. 0. 500000, 0.000000 

1 .QOOOO, 0.000000 

PN: 

0.000000, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

i 

§ 00000, 1.50000. 5.00000 

0.000000, 0.000000, 0.01 0000 

150] 

ee 

I 


a7 


evs 

5 

ee 

if 

XQ: 

on 

Ck: End of SBR comparison 


CE: Comparison with measurements 

CE: Antenna at original position 

CE: containing surrounding walls, floor, ceiling 
CE: ALL exterior windows, overhead beam 
CE: ALL NORTHERN partition walls (v.1, v.2 and v.3) 
CE: with “open” doors (DOUBLE v.3) 

CE AEE SOLITHERN partitions 

CE:ALL the rest of INTERIOR PARTITIONS 
CE:ALL the wooden tables 

CE incudnoine MERARIGBEOGk 

CE: and the interior columns 

CE: using volumetric instead of linear pattern 
UN: units in meters 

l 

US: source units in wavelengths 

() 

FR. trequency of operation 

2.45000 

ZI) CONnGIele Ss PlOpeniies 

-4 

| 

0.061200, 4.00000, 0.00436300, 1.00000, 0.000000 
PM: building's exterior structure 

§ 

0.000000. 0.000000. 0.000000 

0.000000. 0.000000, 3.81000 

0.000000, 39.8272, 3.81000 

0.000000. 39.8272, 0.000000 

39.8272, 39.8272, 0.000000 

50 62 239.8272. )-6 NOOO 

39.8272, 0.000000, 3.81000 

39.8272, 0.000000, 0.000000 

() 


ad io — 
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© GV Vif to oe 


iO. i Cn he Oo, Io i 


di —- OO ~1 HD W 


CE: Building exterior windows(NORTHERN part) 
CE: Exterior Windows 
ZD: glass properties 


| 

0.02000, 4.00000, 0.000600000, 1.00000. 0.000000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 

0.000000. 5.62610, 3.63220 

0.000000. 0.444500, 3.63220 

RR: translating exterior window 
0.000000, 5.68960. 0.000000 

0.000000, 0.000000, 90.0000. 0.000000 


oo 


PG: exterior window 

4,0 

0.000000, 0.444500. 0.177800 
0.000000, 5.62610. 0.177800 
0.000000, 5.62610. 3.63220 
0.000000. 0.444500. 3.63220 
RR: translating extenor window 
0.000000, 5.68960. 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000. 5.62610, 0.177800 
0.000000 55.626 105 365220 

0. VOOOOD0. 0.444500, 3.63220 
RR: translating exterior window 
0.000000. 5.68960, 0.000000 
0.000000. 0.000000. 90.0000, 0.000000 
PG: extenor window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 
0.000000. 5.62610, 3.63220 
0.000000, 0.444500, 3.63220 
RR: translating extenor window 
0.000000. 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: exterior window 

4.0 

0.000000. 0.444500, 0.177800 
0.000000. 5.62610. 0.177800 
0.000000, 5.62610, 3.63220 
0.000000. 0.444500, 3.63220 
RR: translating exterior window 
0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: exterior window 

4,0 

0.000000, 0.444500. 0.177800 
0.000000. 5.62610. 0.177800 

0. QOO0000. 5.62610. 3.63220 

(). OOOO00. 0.444500, 3.63220 
RR: translating exterior window 


100 


0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

PG: exterior window 

4.0 

0.000000, 0.444500, 0.177800 

0.000000. 5.62610, 0.177800 

0.000000. 5.62610, 3.63220 

0.000000. 0.444500, 3.63220 

RR: reverting to onginal(GLOBAL)coordinate system 
0.000000. -34.1376, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

Ci: 

CE:END of building NORTHERN windows 
Cry 

CE:STARTING exterior WESTERN windows 
Gry 

RR: rotating exterior window 

0.000000, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, -90.0000 

PG: extenor window 

4.0 

0.000000. 0.444500. 0.177800 

0.000000. 5.62610. 0.177800 


0.000000. 0.444500. 3.63220 

RR: translating exterior window 
0.000000, 5.68960. 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 
0.000000, 5.62610, 3.63220 

0.000000, 0.444500, 3.63220 

RR: translating exterior window 
0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: exterior window 

4.0 

0.000000, 0.444500, 0.177800 
0.000000. 3.62610. 0.177800 
0.000000. 3.62610. 3.63220 

0.000000, 0.444500. 3.63220 


RR: translating exterior window 
0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: extenor window 

ao) 

QO OOOOO0D. 0.444500. 0.177800 
0.000000. 5.62610, 0.177800 
0.000000. 5.62610. 3.63220 

0.000000. 0.444500. 3.63220 

RR: translating exterior window 
0.000000. 5.68960. 0.000000 
0.000000. 0.000000, 90.0000, 0.000000 
PG: exterior window 

4.0 

O:00C000. 0449500. 02177500 
0.000000, 5.62610, 0.177800 

OU: COG00U% Ss 626101316357 2) 

0.000000, 0.444500, 3.63220 

RR: translating exterior window 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000. 90.0000, 0.000000 
PG: extenor window 

4.0 

0.000000. 0.444500. 0.177800 
0.000000. 5.62610. 0.177800 

0) QOU0G0S 35626110" 3763220 

0.000000. 0.444500. 3.63220 

RR: translating exterior window 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PG: extenor window 

4.0 

Q.QQ00000. 0.444500. 0.177800 
0.000000. 5.62610. 0.177800 
Q.QOOO0O. 5.62610. 3.63220 

0.000000, 0.444500, 3.63220 

CE. 

CE:END of building WESTERN exterior windows 
Gr 

CE: 

CE:STARTING SOUTHERN extenor windows 
CE 


RR: rotaung/translating exterior window 


Q.QO00000. 5.68960. 0.000000 
0.000000, 0.000000, 90.0000, 90.0000 
CM: PG: exterior window 

Givi 470 

CM: 0.000000, 0.444500, 0.177800 
CM: 0.000000, 5.62610, 0.177800 
CM: 0.000000, 5.62610, 3.63220 

CE: 0.000000. 0.444500, 3.63220 

RR: rotatng/translating exterior window 
0.000000. 5.68960, 0.000000 
0.000000. 0.000000, 90.0000, 0.0000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
Q.000000, 5.62610, 0.177800 
O-G00000, 3:62610, 3.63220 

(). 000000, 0.444500. 3.63220 

RR: rotaung/translating exterior window 
OQ.Q000000. 5.68960. 0.000000 
0.000000, 0.000000, 90.0000. 0.0000 
PG: exterior window 

4.0 

0.000000, 0.444500, 0.177800 
0.000000. 5.62610, 0.177800 
0.000000. 5.62610, 3.63220 

0.000000, 0.444500, 3.63220 

RR: rotating/translating exterior window 
0.000000. 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4.0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610. 0.177800 
0.000000. 5.62610. 3.63220 

O.QOOOO000, 0.444500. 3.63220 

RR: rotaung/translating exterior window 
0.000000, 5.68960. 0.000000 
0.000000, 0.000000, 90.0000, 0.0000 
CM: PG: exterior window 

CM 4-0 

CM: 0.000000, 0.444500, 0.177800 
CM: 0.000000. 5.62610, 0.177800 
CM: 0.000000. 5.62610, 3.63220 
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CE: 0.000000. 0.444500. 3.63220 

RR: rotating/translating extenior window 
0.000000, 5.68960. 0.000000 

0.000000, 0.000000, 90.0000. 0.0000 
CM: PG: exterior window 

CM: 4, 0 

CM: 0.000000. 0.444500, 0.177800 
CM: 0.000000, 5.62610, 0.177800 

Git 0000000: 3. 626105363220 

CE: 0.000000, 0.444500, 3.63220 

RR: rotauing/translating exterior window 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
CM: PG: exterior window 

CM: 4, 0 

CM: 0.000000, 0.444500, 0.177800 
CM: 0.000000. 5.62610, 0.177800 

CM: 0.000000. 5.62610. 3.63220 

CE: 0.000000. 0.444500, 3.63220 

Ci 

CE:END of SOUTHERN windows 

Cle: 

CE: 

CE:STARTING EASTERN exterior windows 
Cre 

RR: rotating/translaung exterior window 
0.000000. 5.68960. 0.000000 

0.000000. 0.000000, 90.0000, 90.0000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 

0.000000, 5.62610, 3.63220 

0.000000. 0.444500, 3.63220 

RR: rotaung/translating exterior window 
0.000000. 5.68960. 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4.0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 

0.000000. 5.62610, 3.63220 

0.000000. 0.444500. 3.63220 
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RR: rotating/translating extenor window 
0.000000. 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4,0 

0.000000. 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 

0.000000. 5.62610, 3.63220 

0.000000. 0.444500, 3.63220 

RR: rotating/translating exterior window 
0.000000, 5.68960. 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4,0 

0.000000. 0.444500, 0.177800 
0.000000. 5.62610. 0.177800 

0.000000. 5.62610. 3.63220 

0.000000, 0.444500, 3.63220 

RR: rotating/translating exterior window 
0.000000. 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4,0 

0.000000, 0.444500, 0.177800 
0.000000, 5.62610, 0.177800 

OOCCCOCs G2010Ns 63720 

0.000000, 0.444500, 3.63220 

RR: rotating/translating exterior window 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: extenior window 

4.0 

0.000000, 0.444500. 0.177800 
0.000000. 5.62610, 0.177800 

0.000000. 5.62610. 3.63220 

0.000000, 0.444500. 3.63220 

RR: rotating/translating exterior window 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.0000 
PG: exterior window 

4.0 

0.000000, 0.444500. 0.177800 
0.000000, 5.62610, 0.177800 


ders 


VY OOOOOO. 5.62610. 3.63220 

0.000000. 0.444500. 3.63220 

Cir: 

CE:END of building EASTERN windows 
CE: 

CE: reverting to ORIGINAL (GLOBAL) 
CE: coordinate system (rotating) 

RR: 

0.000000, 0.000, 0.000000 

0.000000, 0.000000. 90.0000, -90.000000 
GE teverine to ORIGINAL (GEOBAE) 
CE: coordinate system (translating) 

RR: 

-5.68960, -39.8272, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
Cr: 

CE:END of building ALL exterior windows 
ele 

CE: partutions between rooms 

ZD: plaster properties 

| 

2 

0.0612000, 2.30000, 0.00174500, 1.00000, 0.000000 
0.0100000. 2.30000, 0.00174533, 1.00000, 0.000000 
Or 

4.0 

0.000000. 2.84480, 0.000000 

0.000000. 2.84480, 3.81000 

5.68960. 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480. 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
Gr 

4.0 

Q QOOOOO. 2.84480. 0.000000 

Q.OVOVO0D. 2.84480, 3.8 1000 

5.68960, 2.84480, 3.81000 

5.68960. 2.84480, 0.000000 

RR: 

0.000000, 2.84480. 0.000000 

0 000000. 0.000000, 90.0000, 0.000000 
Gi: PG: 


CM: 4. 0 

CM: 0.000000, 2.84480, 0.000000 

CM: 0.000000, 2.84480, 3.8 1000 

CM: 5.68960, 2.84480, 3.81000 

CE: 5.68960, 2.84480, 0.000000 

ZD: concrete wall containing metal beams 
! 

0.06 12000, 8.24700. 0.00436300, 1.00000, 0.000000 
0.06 12000. 8.24700. 0.00436300. 1.00000, 0.000000 
Roe 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PGs: 

+.0 

0.000000. 2.84480. 0.000000 

0.000000, 2.84480, 3.8 LOOO 

5.68960. 2.84480. 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
CMR G: 

CM: 4.0 

CM: 0.000000. 2.84480, 0.000000 

CM: 0.000000, 2.84480, 3.8 1000 

CM: 5.68960. 2.84480. 3.81000 

CE: 5.08960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PG: 

4.0 

0.000000. 2.84480, 0.000000 

0.000000. 2.84480. 3.8 1000 

5.68960. 2.84480. 3.8 1000 

5.68960, 2.84480, 0.000000 

CE: change material properties. 

ZD: plaster properties 

| 

2 

0.061 2000, 2.30000, 0.00174500, 1.00000. 0.000000 
0.0 100000. 2.30000. 0.00 174533, 1.00000. 0.000000 


1G: 


RR. 

0.000000. 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
CM: PG: 

GM: 4.0 

CM: 0.000000. 2.84480. 0.000000 

CM. 0.000000. 2.84480. 3.81000 

CM: 5.68960. 2.84480. 3.81000 

CM: 5.68960. 2.84480. 0.000000 

CE 

RR: 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
Ge 

4.0 

0.000000. 2.84480, 0.000000 

0.000000. 2.84480. 3.81000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000. 0.000000 
PG: 

a0) 

0.000000. 2.84480. 0.000000 

(). QOQU0U, 2.84480. 3.8 1000 

5.68960. 2.84480. 3.81000 

5.68960. 2.84480. 0.000000 

Rene 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
RG. 

4.0 

0.000000. 2.84480, 0.000000 
0.000000. 2.84480, 3.81000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

Rue 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PGs 

a8, 

Q. QOUOVO0. 2.84480, 0.000000 


() OOOO0O. 2.84480. 3.8 1000 
5.68960, 2.84480, 3.8 1000 
5.68960, 2.84480, 0.000000 


RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

Eo 

+.0 

0.000000. 2.84480, 0.000000 

0.000000, 2.84480, 3.8 1000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480. 0.000000 

0.000000. 0.000000. 90.0000. 0.000000 

Fic 

+.0 

0.000000, 2.84480, 0.000000 

0.000000, 2.84480, 3.81000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: reverting to global axis system 

0.000000. -34. 1376, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

ZD: concrete wall containing metal beams 

| 

2 

0.06 12000, 5.84000, 0.00141600, 1.00000, 0.000000 
0.06 12000, 5.84000, 0.00141600, 1.00000, 0.000000 
CE: exposed concrete containing beams over 

CE: the interior partitions 

PG: 

4.0 

5.68960. 0.000000, 3.8 1000 

5.68960. 0.000000, 3.40360 

5.68960, 39.8272, 3.40360 

5.68960, 39.8272, 3.81000 

RR: moving the fake partition to the right position 
5.68960, 0, 0.000000 

0.000000. 0.000000, 90.0000, 0.000000 

CM: 

CM: Building the simplified partition (doors open) 


IG ORS, 


CEs 

ZD. plaster properties (prevailing) 

| 

2 

0.0612000, 2.30000, 0.00174500. 1.00000, 0.000000 
0.0100000, 2.30000, 0.00174533, 1.00000, 0.000000 
PG: simplified partition (doors open) 

8. 0 

0.OQ0000. 0.000000, 2.13360 
0.000000. 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000, 4.46720, 2.13360 

0.000000, 5.68960, 2.13360 

0.000000. 5.68960, 3.40360 

0.000000. 0.000000. 3.40360 

RR: 

0.050000. 0.000000, 0.000000 
0.000000. 0.000000. 90.0000, 0.000000 
PG: simplified partition (doors open) 

8, 0 

0.000000, 0.000000, 2.13360 
0.000000, 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 
0.000000. 4.46720, 0.000000 
0.000000. 4.46720, 2.13360 

0.000000. 5.68960, 2.13360 

0.000000, 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

RR: 

-0.050000, 5.6896, 0.000000 

0.000000. 0.000000, 90.0000, 0.000000 
GE 

CE: simplified partition v.3 
CE:(without door, windows etc) 

CEs 

CE: we insert a “double-layer” of 

CE: plaster wall, each layer with thickness 6cm 
CE 

CEs 

ap 

ZD: plaster properties (prevailing) double-layer 
| 


Hee 


5 


0.061 2000, 2.30000, 0.00174500, 1.00000, 0.000000 
0.0612000, 2.30000, 0.00174500, 1.00000, 0.000000 
PG: simplified partition (door open) 

6. 0 

0.000000, 0.000000, 0.000000 

0.000000. 4.46720, 0.000000 

0.000000. 4.46720, 2.13360 

0.000000. 5.68960, 2.13360 

0.000000, 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

GE. 

CE: we now revert to single-layer 

CE. 

ZD: plaster properties (prevailing) 

| 

| 

0.061200, 2.30000, 0.00174533, 1.00000, 0.000000 
CM: RR: 

CM: 0.020000, 0.000000, 0.000000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 
CM: PG: simplified partition (door open) 
CM: 6, 0 

CM: 0.000000, 0.000000, 0.000000 

CM: 0.000000. 4.46720, 0.000000 

CM: 0.000000. 4.46720. 2.13360 

CM: 0.000000, 5.68960. 2.13360 

CM: 0.000000, 5.68960, 3.40360 

CM: 0.000000, 0.000000, 3.40360 

CM: RR: reverting to original axis-system 
CM: -0.020000, 0.000000, 0.000000 

CE: 0.000000, 0.000000, 90.0000, 0.000000 
Cr: 

CE: building the rest of NORTHERN partitions 
Gr: 

RR: moving partition v.1 to the nght 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

RR: moving partition v.1 to the right 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

RR: moving partition v.] to the nght 
0.000000. 5.68960, 0.000000 


0 VOOVOD, 0.000000, 90.0000, 0.000000 
PG: simplified partition (doors open) 
8. 0 

0.000000, 0.000000, 2.13360 
0.000000, 1.22240, 2.13360 
0.000000. 1.22240, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000. 4.46720, 2.13360 
0.000000. 5.68960. 2.13360 
0.000000, 5.68960, 3.40360 
0.000000, 0.000000, 3.40360 

RR: moving partition v.1 to the nght 
0.010000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: simplified partition (doors open) 
8. 0 

0.000000. 0.000000, 2.13360 
0.000000. 1.22240, 2.13360 
0.000000. 1.22240, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000, 4.46720, 2.13360 
0.000000, 5.68960, 2.13360 
0.000000, 5.68960, 3.40360 
0.000000, 0.000000, 3.40360 

RR: moving partition v.1 to the nght 
-Q.010000. 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: simplified partion (doors open) 
Sa0) 

0.000000, 0.000000, 2.13360 
0.000000, 1.22240, 2.13360 
0.000000, 1.22240, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000. 4.46720. 2.13360 
0.000000. 5.68960. 2.13360 
0.000000. 5.68960, 3.40360 
0.000000. 0.000000. 3.40360 

Oe es: 

CE: End of building partitions v.1 on NORTHERN side 
CE: 

Ce 

CE: Starting partition v.2 

Cre. 


Cle 

RR. moving the avis system to the left 
0.000000, -1 7.0688, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
Gr: 

CE: simplified partition v.2 
CE:(without doors, windows etc) 

PG: simplified partition (doors open) 

6. 0 

0. 000000, 0.000000, 2.13360 
0.000000. 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 
0.000000, 5.68960, 0.000000 
0.000000, 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

CE 

Cre. 

CE:END of building NORTHERN partitions 
ole. 

RR: moving back the antenna 
0.000000, -11.3792, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE. 

Cle: 

CE oraning the SOU) BERN pant 

CE: 

RR: moving the axis System 

28.448, -5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
Cie: 

Cr: 

CE: partitions between rooms 

ZD: plaster properties 

| 

| 

0.612000, 2.30000, 0.00174533, 1.00000, 0.000000 
AG 

+.0 

0.000000, 2.84480, 0.000000 
0.000000, 2.84480, 3.81000 

5.68960. 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 


a 


Q.QOOOOO. 2.84480, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
PG: 

4.0 

0.000000, 2.84480, 0.000000 
0.000000. 2.84480, 3.81000 

5.68960. 2.84480. 3.81000 

5.68960. 2.84480. 0.000000 

RR: 

0.000000. 2.84480, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
GMI: PG: 

CM: 4, 0 

CM: 0.000000, 2.84480, 0.000000 

CM: 0.000000. 2.84480, 3.81000 

CM: 5.68960, 2.84480, 3.81000 

CE: 5.68960, 2.84480. 0.000000 

ZD: concrete wall containing metal beams 
2 

0.0612000, 5.84000, 0.00141600, 1.00000, 0.000000 
0.0612000, 5.84000, 0.00141600, 1.00000, 0.000000 
Re 

0.000000. 2.84480. 0.000000 

() OCOOOOO. 0.000000. 90.0000. 0.000000 
CMI, PG: 

CM: 4.0 

CM: 0.000000, 2.84480, 0.000000 

CM: 0.000000, 2.84480, 3.81000 

CM: 5.68960, 2.84480, 3.81000 

CE: 5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480. 0.000000 
0.000000. 0.000000. 90.0000. 0.000000 
rie 

(0) 

0.000000, 2.84480. 0.0OO0000 
Q.0Q0000, 2.84480, 3.81000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000. 2.84480. 0.000000 

Q.Q00000. 0.000000. 90.0000. 0.000000 


Pon 

4.0 

0.000000, 2.84480, 0.000000 
0.000000. 2.84480, 3.8 1000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

CE: change material properties 

ZD: plaster properties 

| 

| 

0.612000, 2.30000, 0.00174533, 1.00000, 0.000000 
RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
eG: 

CM: 4. 0 

CM: 0.000000. 2.84480. 0.000000 

CM: 0.000000, 2.84480. 3.8 L000 

CM: 5.68960. 2.84480, 3.81000 

CE: 5.68960, 2.84480, 0.000000 

i oh 

RR: 

0.000000, 2.84480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PGs 

4,0 

0.000000. 2.84480, 0.000000 
0.000000, 2.84480, 3.8 1000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480. 0.000000 

0.000000, 0.000000. 90.0000, 0.000000 
PG: 

4.0 

0.000000, 2.84480, 0.000000 
0.000000, 2.84480, 3.8 1000 

5.68960, 2.84480, 3.81000 

5.68960, 2.84480, 0.000000 

Rie 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
PG: 


AS 


4.0 

0.000000. 2.84480. 0.000000 
Q.OOOOOO. 2.84480, 3.8 1000 
5.68960, 2.84480, 3.81000 
5.68960, 2.84480, 0.000000 

REe 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
PG: 

4.0 

0.000000, 2.84480, 0.000000 
0.000000, 2.84480, 3.8 1000 
5.68960, 2.84480, 3.81000 
5.68960, 2.84480, 0.000000 

RR: 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
Gs: 

4.0 

0.000000, 2.84480. 0.000000 
0.000000, 2.84480, 3.81000 
5.68960, 2.84480, 3.8 1000 
5.68960, 2.84480, 0.000000 

Rie 

0.000000, 2.84480, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
P&: 

4,0 

0.000000, 2.84480, 0.000000 
0.000000, 2.84480, 3.8 1000 
5.68960, 2.84480, 3.8 1000 
5.68960, 2.84480, 0.000000 

RR: reverting to global axis system 
-5.6896, -34.1376. 0.000000 
Q.0O00000, 0.000000. 90.0000, 0.000000 
Ci PG: 

CM: 4.0 

CM: 0.000000. 1.5, 0.000000 

CM: 0.000000, 1.5, 3.81000 

CM 5 0s9CUr >. so 00D 

CE: 5.68960, 1.5, 0.000000 

CE 

GE. 


Cr: 

Cle 

CE 

CE: Burlding the partitions v.1,v.2,v.3 

Oe ihe le iticiVN Side 

GE: 

CE 

GE; 

CE: 

GE: 

ZD: concrete wall containing metal beams 

0.06 12000, 5.84000, 0.00141600, 1.00000, 0.000000 
0.0612000, 5.84000, 0.00141600, 1.00000, 0.000000 
CE: exposed concrete containing beams over 

CE: the interior partitions 

Ge 

4.0 

5.68960, 0.000000, 3.8 1000 

5.68960, 0.000000, 3.40360 

5.68960. 39.8272, 3.40360 

5.68960, 39.8272, 3.81000 

RR: moving the fake partition to the right position 
5.68960, 0, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 

CM: 

CM: Building the simplified parution (doors open) 
CE: 

ZD: plaster properties (prevailing) 

| 

= 

0.061200, 2.30000, 0.00174533, 1.00000, 0.000000 
PG: simplified partition (doors open) 

8, 0 

0.000000, 0.000000, 2.13360 

0.000000, 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 

0.000000, 4.46720, 0.000000 

0.000000. 4.46720, 2.13360 

0.000000, 5.68960, 2.13360 

0.000000, 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 


ssa 


Ronse 

0.050000. 0.000000, 0.000000 
0.000000. 0.000000, 90.0000. 0.000000 
PG: simplitied parutuon (doors open) 
$. 0 

0.000000. 0.000000, 2.13360 
0.000000. 1.22240, 2.13360 
0.000000. 1.22240, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000. 4.46720, 2.13360 
0.000000. 5.68960, 2.13360 

0 OOOOOO. 5.68960. 3.40360 
0.000000. 0.000000, 3.40360 

RR: 

-0.050000. 5.6896, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
Cr: 

CE: simplified partition v.3 
CE:(without door, windows etc) 

PG: simplified partiuon (door open) 
6.0 

0.000000, 0.000000. 0.000000 
0.000000, 4.46720, 0.000000 
0.000000, 4.46720, 2.13360 
0.000000. 5.68960, 2.13360 
0.000000, 5.68960. 3.40360 
0.000000, 0.000000, 3.40360 

RR: 

().020000, 0.000000. 0.000000 
0.000000, 0.000000, 90.0000. 0.000000 
PG: simplified partion (door open) 
6,0 

0.000000, 0.000000, 0.000000 
0.000000, 4.46720, 0.000000 
0.000000, 4.46720. 2.13360 
0.000000, 5.68960. 2.13360 

Q). OOOO. 5.68960. 3.40360 

0 OOOQO0O, 0.000000. 3.40360 

RR: reverting to orginal axis-system 
-0.020000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE. 


138 


CE: building the rest of NORTHERN partitions 
Cs 

RR: moving partition v.] to the nght 
0.000000. 5.68960. 0.000000 
0.000000. 0.000000. 90.0000. 0.000000 
RR: moving partition v.1 to the nght 
0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
RR: moving partition v.1 to the nght 
0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: simplified partition (doors open) 

8. 0 

0.000000, 0.000000, 2.13360 

0.000000, 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 

0.000000, 4.46720, 0.000000 

0.000000, 4.46720, 2.13360 

0.000000, 5.68960, 2.13360 

0.000000. 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

RR: moving parttuon v.{ to the nght 
0.000000. 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PG: simplified partition (doors open) 
520 

0.000000, 0.000000, 2.13360 
0.000000, 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 

0 .QOOQ000, 4.46720. 0.000000 

0.000000. 4.46720, 2.13360 

0.000000. 5.68960, 2.13360 

0.000000, 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

RR: moving partition v.1 to the nght 
0.000000, 5.68960, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PG: simplified partition (doors open) 

§. 0 

0.000000, 0.000000, 2.13360 

0.000000, 1.22240, 2.13360 

0.000000, 1.22240, 0.000000 
0.000000, 4.46720. 0.000000 


}-- 


f 3 


LO 


0.000000, 4.46720, 2.13360 

0.000000, 5.68960, 2.13360 

0.000000. 5.68960. 3.40360 

0.000000. 0.000000, 3.40360 

is. 

CE: End of building partitions v.1 on NORTHERN side 
CE 


Cle 
CE: Starting partition v.2 
Cr 
CE; 


RR: moving the axis system to the left 
0.000000. -17.0688. 0.000000 
0.000000. 0.000000. 90.0000, 0.000000 
GE 

CE: simplified partition v.2 
CE:(without doors, windows etc) 

PG: simplified partition (doors open) 
6.0 

0.000000, 0.000000, 2.13360 
0.000000, 1.22240, 2.13360 

0.000000, |.22240, 0.000000 
0.000000, 5.68960, 0.000000 
0.000000. 5.68960, 3.40360 

0.000000, 0.000000, 3.40360 

CE 

Cds) 

CE. 

Oley 

GE. 

RR: reverting to GLOBAL coordinate system 
-34.1376, -17.0688, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE 

CE 

CE 

CE: Building the rest of interior partitions 
CE 

Cr. 

CE 

CE 

RR: translating the four NORTH-SOUTH partitions 
5.68960, 11.3792, 0.000000 


0.000000. 0.000000. 90.0000, 0.000000 
PG: interior partition (NORTH-SOUTH) 
§. 0 

0.000000, 0.000000, 0.000000 
0.000000, 0.000000, 3.8 1000 

3.04800, 0.000000, 3.81000 

3.04800, 0.000000, 0.000000 

2.43840, 0.000000, 0.000000 

2.43840, 0.000000, 2.13360 

0.914400. 0.000000, 2.13360 
0.914400. 0.000000. 0.000000 

CE: translating the second out of the four 
CE: NORTH-SOUTH partitions 

RR: 

0.000000, 5.68960, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: intenor partition (NORTH-SOUTH) 
8, 0 

0.000000, 0.000000. 0.000000 
0.000000. 0.000000. 3.8 1000 

3.04800, 0.000000, 3.81000 

3.04800. 0.000000, 0.000000 

2.43840, 0.000000, 0.000000 

2.43840, 0.000000, 2.13360 

0.914400, 0.000000, 2.13360 
0.914400, 0.000000, 0.000000 

CE: translating the third out of the four 
CE:NORTH-SOUTH partitions 

RR: 

250.0) UOOG00 

0.000000. 0.000000, 90.0000, 0.000000 
PG: interior partition (NORTH-SOUTH) 
8, 0 

0.000000, 0.000000, 0.000000 
0.000000, 0.000000, 3.8 1000 

3.04800, 0.000000. 3.8 1000 

3.04800, 0.000000, 0.000000 

2.43840. 0.000000. 0.000000 

2.43840. 0.000000, 2.13360 

0.914400. 0.000000. 2.13360 
0.914400, 0.000000, 0.000000 

CE: translating the forth out of the four 
CE: NORTH-SOUTH partitions 


ib il 


ine 

0. -5.6896, 0.000000 

0.000000. 0.000000, 90.0000, 0.000000 
PG: intenor purution (NORTH-SOUTH) 
S. 0 

0 000000. 0.000000. 0.000000 
0.000000, 0.000000. 3.8 1000 

3.04800, 0.000000. 3.81000 

3.04800, 0.000000, 0.000000 

2.43840, 0.000000, 0.000000 

2.43840, 0.000000, 2.13360 

0.914400. 0.000000, 2.13360 
0.914400. 0.000000, 0.000000 

CE 

CE 

CE 

RR: reverting to GLOBAL coordinates 
-31.0796, -11.3792, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE 

CE 

CE: Building the WEST-EAST partition 
CE. 

RR: translaung the WEST-EAST partition 
17.0688, 0.000000. 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PG: WEST-EAST interior partition 

8. 0 

0.000000, 0.000000, 0.000000 
0.000000. 0.000000, 3.8 1000 

0.000000. 11.3792, 3.81000 

0.000000, 11.3792, 0.000000 
0.000000, 10.6172, 0.000000 
0.000000, 10.6172. 2.13360 

0.000000. 9.09320, 2.13360 

0.000000, 9.09320, 0.000000 

RR: reverung to GLOBAL coordinates 
-17.0688. 0, 0.000000 

0.000000. 0.000000. 90.0000. 0.000000 
os 

Car, 

iE: 

CE: Building central room’s light long 


Ie2 


CE: NORTH-SOUTH partition 

CE. 

RR: translating the LONG central partition 
11.3792, 11.3792, 0.000000 

0.000000, 0.000000, 90.0000. 0.000000 
PG: central room’s long partition 

ee 

0.000000, 0.000000. 0.000000 

0.QO00000. 0.000000. 3.8 1000 

17.0688, 0.000000, 3.8 1000 

17.0688. 0.000000, 0.000000 

14.2240, 0.000000, 0.000000 

14.2240, 0.000000, 2.13360 

11.3792, 0.000000, 2.13360 

11.3792, 0.000000, 0.000000 

5.68960, 0.000000, 0.000000 

5.68960, 0.000000, 2.13360 

+.77520, 0.000000, 2.13360 

4.77520, 0.000000, 0.000000 

CE: 

Cr, 

CE: 

RR: reverting to GLOBAL coordinates 

se ke oes es eee a erm one 010.010.8 

0.000000, 0.000000, 90.0000, 0.000000 
Cie 

Cie: 

Cr 

CE: Building the concrete wall in front of 
CE: the stairs; using “artificial” 

CE: constitutive parameters 

ZD: concrete wall containing metal beams 
0.06 12000, 8.247, 0.00141600, 1.00000, 0.000000 
0.06 12000, 8.247, 0.00141600, 1.00000, 0.000000 
RR: translating the concrete wall 
8.73760, 11.3792, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
CE: concrete wall in front of the stairs 
PM: 

16 

0.00500000. 0.00500000, 0.000000 


eS 


J VOSOOQVO0. O.0OOS500000. 3.8 1000 
0. OOSO0000. 5.68960, 3.81000 
0.00500000. 5.68960, 0.000000 
2.64160, 5.68960, 0.000000 
2.64160, 5.68960, 3.81000 
2.64160, 5.65960, 3.81000 
2.64160. 5.65960, 0.000000 
0.0300000, 5.65960, 0.000000 
0.0300000. 5.65960. 3.8 1000 
0.0300000, 0.0300000. 3.8 1000 
0.0300000, 0.0300000. 0.000000 
2.64160, 0.0300000, 0.000000 
2.64160, 0.0300000, 3.81000 
2.64160, 0.00500000, 3.81000 
2.64160, 0.00S500000, 0.000000 
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CE: 

RR: reverting to GLOBAL coordinates 
-§.73760, -11.3792, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
or: 

GE. 

Ge 

Rk translanne ihe Vile Sit 

11.4792. 17.0688. 0.000000 

0.000000, 0.000000. 90.0000. 0.000000 
CE 

CE 

CE: MESH of the instrument room 

GE: 

CE: 

ZD: using made-up properties for MESH 
| 

| 

0.0100000,. 5.896, 0, 1.00000, 0.000000 
PG: metallic MESH substitute 

4,0 

0.000000, 0.000000, 0.000000 
0.000000, 0.000000, 3.8 1000 

16.9088. 0.000000. 3.8 1000 

16.9088. 0.000000. 0.000000 

RR: reverting to GLOBAL coordinates 
-}}].4792, -17.0688. 0.000000 


0.000000. 0.000000, 90.0000, 0.000000 

CE 

CEs 

O) 2 

CE: changing properties again 

CE: 

Cie 

ZD: concrete wall containing metal beams 

0.061 2000, 8.247, 0.00141600, 1.00000, 0.000000 
0.0612000, 8.247, 0.00141600, 1.00000, 0.000000 
RR: translating the second and last concrete wall 
28.4480, 11.3792, 0.000000 

0.000000. 0.000000, 90.0000. 0.000000 

CE: second and last concrete wall 

Gre. 

CE: 

PM: 

8 

0.00500000, 0.00500000, 0.000000 
0.00500000, 0.00500000, 3.8 1000 

0.00500000, 5.68960, 3.81000 

0.00500000. 5.68960, 0.000000 

2.64160. 5.68960, 0.000000 

2.64160. 5.68960. 3.81000 

2.64060, 0.00500000, 3.81000 

2.64060. 0.00500000, 0.000000 
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SS 


RR: reverting to GLOBAL coordinates 
-28.4480. -11.3792, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE, 

Cr. 

CE: 

a 

ee: 

CE: Building the wooden tables 


ZD: wood properties 

l 

| 

0.06 12000, 1.50000, 0.0100000, 1.00000, 0.000000 
GE: 

Ce 

RR: translating the tables v.1 

7.82320, 0.914400. 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PG: wooden table v.1 

4.0 

0.000000, 0.000000, 0.914400 
0.000000, 8.22960, 0.914400 

1.82880, 8.22960. 0.914400 

|.82880. 0.000000. 0.914400 

RR: translating the second table v. 1] 

2.5. 0. 0.000000 

0.000000, 0.000000. 90.0000, 0.000000 
PG: wooden table v.1] 

4,0 

0.000000, 0.000000, 0.914400 
0.000000. 8.22960. 0.914400 

1.82880. 8.22960, 0.914400 

1.82880, 0.000000, 0.914400 


a 


RR translating the third table v.1 

2.5. 0. 0.000000 

0.000000. 0.000000. 90.0000. 0.000000 
PG: wooden table v.1 

4,0 

0.000000, 0.000000, 0.914400 
0.000000, §.22960, 0.914400 

| .§2880. 8.22960, 0.914400 

1. §2880, 0.000000, 0.914400 

Oe) 

CE 

RR: reverting to GLOBAL 

-12.82320, -0.914400, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
Oe: 

CE: 

CE: 

CE: Building the wooden tables v.2 
CE: 

GE 

CE. 

RR: translating the wooden tables v.2 
29.6164, 22.7584, 0.000000 

0.000000. 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (Ist of Ist row) 
S 

0.000000. 0.000000, 0.990600 
0.000000. 5.28320, 0.990600 

ea 20 255200 000 

1.47320, 0.000000, 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000, 1.16840 
0.000000, 5.28320, 1.16840 

1.47320, 5.28320, 1.16840 

0 
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iE 

RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (2nd of Ist row) 
§ 

0.000000. 0.000000. 0.990600 
0.000000, 5.28320, 0.990600 

1.47320, 5.28320, 0.990600 

1.47320, 0.000000, 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000, 1.16840 
0.000000, 5.28320, 1.16840 

lea 7520. 5.25520. il6s40 

0 


op ae” a re ae SS ee 


R: translating tables v.2 

-3. 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (3rd of Ist row) 
8 

0.000000, 0.000000, 0.990600 
0.000000, 5.28320, 0.990600 

1.47320. 5.28320, 0.990600 

1.47320. 0.000000, 0.990600 

1.47320. 0.000000. 1.16840 

0 QOOOOO. 0.000000, |. 16840 
0.000000. 5.28320, 1.16840 

Pa 32095. 250207 1 ose 
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0 

RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (4th of Ist row) 
8 

0.000000, 0.000000. 0.990600 
0.000000, 5.28320, 0.990600 

1.47320. 5.28320, 0.990600 

1.47320. 0.000000, 0.990600 

1.47320. 0.000000, 1.16840 

0.000000. 0.000000, 1.16840 
0.000000, 5.28320, 1.16840 

1.47320, 5.28320, 1.16840 

0 
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AoOocwn Oa ON © OU So a =] S12 a ay Se eee 


R: translating tables v.2 

-3, 0.000000. 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (Sth of Ist row) 
S 

0.000000. 0.000000, 0.990600 
0.000000, 5.28320, 0.990600 

147320, 5, 28520,,0, 990600 

1.47320, 0.000000, 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000, 1.16840 
0.000000. 5.28320, 1.16840 

1.47320. 5.28320. 1.16840 
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RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (6th of Ist row) 
§ 

0.000000. 0.000000. 0.990600 
0.000000. 5.28320, 0.990600 

pa 520. 3.283200 990600 

1.47320. 0.000000. 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000, 1.16840 


0.000000. 5.28320. 1.16840 
1.47320, 5.28320. 1.16840 
O 
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CE: 

CE: End of Fikes ROW ol tables v2 

CE: 

RR: translating the SECOND ROW of the tables v.2 
15.0000, 10.3124, 0.000000 

0.000000. 0.000000. 90.0000. 0.000000 

PM: wooden tables v.2 (lst of 2nd row) 

§ 

0.000000, 0.000000, 0.990600 


0.000000, 5.28320, 0.990600 
Leie20, 5.28320, 0.990600 

| 47320, 0.000000, 0.990600 
1.47320. 0.000000, 1.16840 
0.000000, 0.000000, 1.16840 
0.000000, 5.28320, 1.16840 
1.47320, 5.28320, 1.16840 
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CE: 

RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 


PM: wooden tables v.2 (2nd of 2nd row) 
S 

0.000000. 0.000000. 0.990600 
0.000000. 5.28320. 0.990600 
1.47320. 5.28320, 0.990600 
1.47320. 0.000000. 0.990600 
1.47320, 0.000000, 1.16840 
0.000000, 0.000000, 1.16840 
0.000000, 5.28320, 1.16840 
1.47320, 5.28320, 1.16840 

0 
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0 
RR: translating tables v.2 


ibe Ze 


-3, 0.000000, 0.000000 

0.000000, 0.000000. 90.0000, 0.000000 
PM: wooden tables v.2 (3rd of 2nd row) 
§ 

O.Q000000. 0.000000, 0.990600 
0.000000. 5.28320, 0.990600 

1.47320. 5.28320. 0.990600 

1.47320, 0.000000, 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000, 1.16840 

0.000000, 5.28320, 1.16840 

4g 320552652051 o64e 
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RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000. 0.000000. 90.0000, 0.000000 
PN1. wooden tables v.2 (4th of 2nd row) 
S 

0.000000, 0.000000, 0.990600 
0.000000, 5.28320, 0.990600 

1.47320, 5.28320, 0.990600 

1.47320, 0.000000, 0.990600 

1.47320. 0.000000. !.16840 

0.000000, 0.000000. 1.16840 

0.000000. 5.28320. 1.16840 

1.47320. 5.28320. 1.16840 

0 
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RR: translating tables v.2 

-3, 0.000000, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
PM: wooden tables v.2 (Sth of 2nd row) 
§ 

0.000000, 0.000000, 0.990600 
0.000000, 5.28320, 0.990600 

1.47320, 5.28320, 0.990600 

1.47320, 0.000000, 0.990600 

1.47320, 0.000000, 1.16840 

0.000000, 0.000000. 1.16840 

0.000000. 5.28320. 1.16840 
}47320-9-28520. 116540 

0 
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R: translating tables v.2 

-3, 0.000000. 0.000000 

Q.O00000, 0.000000. 90.0000. 0.000000 
PM: wooden tables v.2 (6th of 2nd row) 
§ 

0.000000, 0.000000, 0.990600 
0.000000, 5.28320, 0.990600 

1.47320, 5.28320, 0.990600 

1.47320. 0.000000, 0.990600 

1.47320. 0.000000, 1.16840 

0.000000. 0.000000. 1.16840 

0.000000. 3.28320. 1.16840 

1.47320, 5.28320, 1.16840 

0 
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() 

Cie 

CE:END OF SECOND ROW 

Cie: 

RR: reverting to GLOBAL coordinates 
-14.6164, -33.0708, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
CE: 

GE: 

ClE\ yD OF ALE WOOBEN TABIEES 
Cr: 

Cr: 

CE 

CE: 

CE: Building the metallic block 

Cie 

Giz: 

RR: translating the metallic block 
22.7584, 28.4480, 0.000000 

0.000000, 0.000000, 90.0000, 0.000000 
ZD: metallic properties (PEC) 

0 

PM: Metallic Block 

8 

-13.4112. 4.26720. 0.000000 
(.000000. 4.26720. 0.000000 
0.000000. 0.000000. 0.000000 
-13.4112, 0.000000. 0.000000 
le ao O LO ee 

0.000000, 4.26720. 1.98120 

0.000000, 0.000000, 1.98120 
-13.4112, 0.000000, 1.98120 

0 
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0 

CE. 

we 

Ole 

RR: reverting to GLOBAL 

-22.7584. -28.4480. 0.000000 
0.000000. 0.000000, 90.0000. 0.000000 
Cre: 

GE. 

Gls: 

CE: 

CE: 

ZD: concrete columns containing metal beams 
Es) 

| 

0.061 2000, 4.00000, 0.00436300, 1.00000, 0.000000 
CE: Building the columns 

CE 

CE, 

RR: translating the columns 

go 2 2 WO OU0000 

0.000000, 0.000000, 90.0000, 0.000000 


Cry 

CE: columns 
CE 

PM lirst column 
S 


142 


0.600000, 0.600000, 0.000000 
-0.600000, 0.600000, 0.000000 
-0.600000, -0.600000, 0.000000 
0.600000, -0.600000, 0.000000 
0.600000, 0.600000, 3.8 1000 
-0.600000. 0.600000. 3.8 1000 
-0.600000. -0.600000, 3.8 1000 
0.600000. -0.600000. 3.8 1000 

QO 
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0 

RR: translating the 2nd column 
5.68960, 0.000000, 0.000000 
0.000000, 0.000000. 90.0000, 0.000000 
PM: second column 

8 

0.600000, 0.600000, 0.000000 
-0.600000, 0.600000, 0.000000 
-Q.600000, -0.600000, 0.000000 
0.600000, -0.600000, 0.000000 
0.600000, 0.600000, 3.81000 
-Q.600000, 0.600000, 3.81000 
-Q.600000, -0.600000, 3.8 1000 
0.600000. -0.600000. 3.8 1000 
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RR: translating the 3rd column 
6.5. 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
PM: third column 

S 

0. GO0000. 0.600000. 0.000000 
-Q.600000, 0.600000, 0.000000 
-0.600000, -0.600000, 0.000000 
0.600000, -0.600000, 0.000000 
0.600000, 0.600000, 3.81000 
-0.600000, 0.600000, 3.8 L000 
-Q.600000, -0.600000. 3.8 1000 
0.600000, -0.600000, 3.8 1000 
0 
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RR: translating the 4th column 
4.8792. 0.000000. 0.000000 
().000000. 0.000000. 90.0000. 0.000000 
PM: fourth column 

§ 

0.600000, 0.600000, 0.000000 
-Q.600000, 0.600000, 0.000000 
-0.600000, -0.600000, 0.000000 
0.600000, -0.600000, 0.000000 
0.600000, 0.600000, 3.8 1000 
-0.600000. 0.600000, 3.81000 
-Q.600000. -0.600000, 3.81000 
0.600000, -0.600000, 3.8 L000 
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~_— 
— 


Ly hs 


Oo + mM 


Q 

RR: reverting to GLOBAL coordinates again 
-2§.448. -29.1. 0.000000 

Q QOOOOOO. 0.000000. 90.0000. 0.000000 

CE: 

Oy 

CE:END of building the columns 

Oe: 

Glee 

CE: 

CESEND OF BUILDING THE WHOLE ELOGK 
CE: 

CE. 

Cir: 

SG: 

3.00000, 8.00000, 1 .QQ0000 

0.000000, 0.000000, 90.0000, 0.000000 

-], 0.500000, 0.000000 

1. OOOOO, 0.000000 

CE: For location #2 

CE: uncomment selection out the 

CE: SG command below and 

CE: comment selection out the above 

GCE SG command 

CM: SG: 

CM: 19.9136, 19.9136, 1.00000 

CM: 0.000000, 0.000000, 90.0000, 0.000000 
CM: -1, 0.500000, 0.000000 

CE: |.00000. 0.000000 

CE: For location #3 

CE uncomment seleciicn outline 

CE: SG command below and 

CE: comment selection out the above 

CM: SG: 

CM: 3.00000, 35.5000, | .QOOO00 

CM: 0.000000. 0.000000. 90.0000. 0.000000 
CM: -1. 0.500000. 0.000000 

CE 1.00000, 0.000000 


}- 
i> 
OV 


CE: For location #4 

CEsuncomment selection out the 

CE: SG command below and 

CE: comment selection out the above 
Cries G- 

CM: 19.9136, 38.0000, 1.00000 

CM: 0.000000. 0.000000, 90.0000, 0.000000 
CM: -!. 0.500000. 0.000000 

CE: |.00000, 0.000000 

CE: For location #5 

CE: uncomment selection out the 

CE: SG command below and 

CE: comment selection out the above 
Cri sSG: 

CM: 37.5000, 4.00000, 1.00000 

CM: 0.000000. 0.000000. 90.0000, 0.000000 
CM: -1. 0.500000. 0.000000 

CE: 1.00000, 0.000000 

CE: 

VN: volumetric pattern 

0.000000, 0.000000, 0.000000 
0.000000, 0.000000, 90.0000, 0.000000 
a 

0.244800, 0.244800, 1.50000 
0.244800. 0.244800, 0.000000 

alo? 162 

CE: The following commented-out 

CE: volumetric pattern was used in the 
CE: comparison input file and the 

CE: rest of the input file remained the same 
CM: VN: volumetric pattern 

CM: 0.000000, 0.000000, 0.000000 
CM: 0.000000, 0.000000. 90.0000. 0.000000 
Cle JE 

CM: 1.32500. 1.32000, 1.50000 

CM: 1.39000, 1.24000, 0.000000 

GE sl, 2429 

BIE USE 

CE mt 

Ele 

I 

TF ea 

le 


nO. 
EN: 
CE: End of input file used in the comparison with measurements 
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APPENDIX B. MATHCAD FILES 


In what follows, all MathCAD files used in the thesis calculations and plots generation 
are given in the order they were presented in the main body. 


MathCAD file used in the cylinder room example. 
Reading the data: 

EE = READPRN Cylinder Koom?_ Ex) 

Calculating the total E-field: 





EA. dB C= EE! 
(5) 

EV dB ‘= EE 
(7) 

E,dB ‘= EE 

es 
EX dB 
> 
E, := 10 a 





i= eee 


The Cylinder’s center of mass 1s at distance 7 meters. 


Room with cylinder 


total | tiedd (dB V/in) 





0 2 4 


See en 


6 8 10 liz al 
distance in meters 


End of cylinder room example 


MathCAD file used in the typical room example (comparison) 
Read the Data 


room_data:= READPRN("Marhefka_room?_E.txt)' 
itee data = READERIN Marheikasroome eircewr- xt 
M := rows(room_data) 

M = 1.728x 10° 


We = Cols( room: @ata 


N= 9 
“stan 
2 
ee 10 ey 
ooo 
os 
room. dati. 
9 
oe ask) ae 
oe 
eer 
20 
Ey := 10 


eo), 


—-- -——. 
(3) 


; a 
{ree data 
20 


free_data 5) 


20 


I 
=> 


en: 
(7) 


free_data 
9 
eo (= LO a 


Calculating the total E-field in the room and in free space 





Each step is 6.55715 inches (16.6551cm) long 
N.x 1s the # points in the outer loop 


(Etot) = 1.728 x 10° 


=. 
WY % Fane: 
ING s tS evoten A FA 


TOWS 












; mh 
Ss 


Ny = 32 
q:=0..Ny—1 
Starlg ‘= q:Nx 


StOPg ‘= Starlg + Nx — | 


_ <q 
room 


‘= submatrit Exo, , Slat lg , StOpg, 0, 0} 
a 
Eqo1_dB *= 20-log{ Exoi) 


room_dB = submatrix Ero dB , Slaltg , StOPg .0, 0} 


Max fOOmedis))—— son 3 

mim room_dB) = 0.797 

max Heldudb:—= i maxtiGOnieds) 

Accounting for the 30dB dynamic range of NEC-BSC 
dyn_range:= 30 


AE ol 


tield_limi= max_field_dB- dyn_range 
p := 0.. rows(room_dB) — 1 
room_tield g := if(room_dBp q > field_limitoom_dBp g,field_lim}| 


cous = submatriy Ejo, . Slag. StOPg . 9. 0) 


freed = subImMaia {etot . Slarlg , SLOPg .0. 0) 
=» 

Crot_dB -= 70-log(eto) 

free_dB? = submatriXeio._dB- slarlg , SLOpg 0, 0) 


Max ineeadisja— 56.04 

mimiree -dis)je— 20.55 

max_field_dB= max(free_dB) 

dyn_range:= 30 

field_limit= max_field_dB- dyn_range 

held _limeé 28.04 

p := Q.. rows(free_dB) — | 

free_field q := if(free_dBy gq > field_limitfree_dBp q, field_lim)t 

free iN? = submatrix eyo, Stary , Stopg 0,0] 

dilf iree_room= free ine room: lit 

diff_tree_room_dB q = if( diff_free_roomp |g > 0. 20-log{ diff_free_roorp, q) 0) 
Mile val" Mita tatree ToGml cds 

min_val= —27.444 

diff_fr_dB gq = if( diff_free_room_dB = 0, min_valdiff_free_room_dB , g) 





room_ field 


3 





a 2U 10 0 


room_field 
End of typical room example (comparison) 


MathCAD file used in the typical room example (random variable calculation) 


Calculating the normalized difference random variable in the typical room example 


IES, 


In this program we are aiming at calculating all the random variables of the normalized 
difference between the rooms and the free space field strength (Etot and Etot_free 
respectively). 

We start by reading the associated data files generated by NEC_BSC. 

NEC_BSC_data:= READPRN"Marhefka_room2_E.txt" ) 

M := rows( NEC_BSC_data) 


Nea 1.728 < 10° 


N= COINGNE Gers scala) 

N= 9 

oor = NEC_BSC_data® 

yee NEC_BSC_data'!” 
(2) 


Zeoord = NEC_BSC_data 





; _ io : ’ ? 9 
(coord 7 4.572) a (Scourd — 5.486)" + wares Sale 19) 
free_space_data := READPRN “Marhetka_room?_tree_E.txt" ) 


MM) .= rows(tree_space _data ) 


NIM = 1.728 10° 

NN := cols(free_space_data ) 

NN =9 

rn 

. NEC_BSC_data® 
%” 

10 ee 


\ EC_BSC data’ 


NEC_BSC_ data” 


20 





rows( Ex, ) = 1.728% 10" 





(3) 5 
| tree_space_data 


| 20 
Ex_free ‘= | 10 


(5) 


free_space_data 


20 
Ey tree = 10 
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a =" 
tree_space_data 


20 
EZ tree ‘= 10 


a ae 
ELot_free = (Ecce us (Ercan a (See 
Let us now create the two datamatrices, namely the "datamatrix" which consists of the 
room's field strength values for every observed point and the "freespacedataOl" which 
includes the corresponding data for the free space case. 


Cat air tx = augment (p Leen) 





treespacedataQ! := augment(p Reece) 

At this point, owing to the non-uniform layout of the distance (from the antenna) points in 
the array, we will have to distinguish them into several different groups using the 
following approach: 

lower := floor( min(p )) 

upper := ceil(max{p )) 


lever — 0 

upper = 7 

bins := upper-2 

bins = 14 

—_ upper — lower 
bins 

n= 05 

n :=Q.. bins 

kh -=lower + a-h 


i 


it should be noted that the number of bins may be chosen 1n a fashion which best suits our 
desirable "distance steps’. For instance, should we choose to have groups points separated 
by 50cm, where the maximum distance from the antenna is 49m, our intuitive selection 
should be 98. 

We next introduce a new function, which shall be utilized, several times in our 
subsequent calculations. As input, it expects a two-column matrix and an array. It looks 
for those values of the first column (i.e., in our problem, the distance array 7 7 which fall 
within our made-up groups and then creates a new matrix that contains all those second 
column’s values whose respective first column value satisfy our requirement. It finally 
returms this new matnx which, it should be underlined, normally has a different number 
of elements in each row. 


Niger, 


peIRON bmivecior jaanl a ceane 
Be me ‘ 
k — binvector 
tor |e O.. rows(k) - 2 

D, — 0 
for 1€ 0.. rows(k) - 2 

for }¢€ 0..rows(M) - 1 

D, — augment(B..D,) if K re a < Ky 

tor 1€ 0.. rows(D) — | 

D — submatrix{ D, .0.rows{D } - 1.0. cols(D)] - 4 if D. =i) 


D 
It is high time we used this helping function. 
datamatrix02:= build(datamatrix. k) 
treespacedata02 := build(treespacedataO! .k) 
The following new function is, seemingly, of high importance since it is the one which 
calculates the "normalized difference” random variables. It receives two matrices of the 
sume dimensions and retums 4 new matrix. which in each row contains the values of the 


random vanable pertinent to that rows “distance group.” 


aimenence( x. 1) = Vitor ke O.rows( xX) —1 
C. — 0 


i for je 0..rows( X) - I 


iets cols( X ] 


A Xx 


i} 
Be Y. 
J 
ic 0 
for te Van 
= Ag i ~ Bo 
ditt. — 
Bo 


ne cols( X 
c — stack(C,, dift} 
©, — (c)" 
ke rows(C.) 
— cols(C ] 
pC submairn{C 00K - 111 1} 


C 
We now get rid of the redundant first row zeros and then calculate the normalized 
difference variables. 
new_datamatri1x02:= submatrix(datamatrixO2, 1, rows(datamatrix02) - 2,0, cols(datamatrix0O2 —- 1) 
new_freespacedata02 := submatrix(freespacedata02 , 1, rows(freespacedata02 ) — 2,0, cols(freespacedataQ2 ) — 1) 
norm_diff := difference(new_datamatrix02, new_freespacedataQ2 ) 
Surprisingly and unexpectedly the Mathcad appears to suffer form a bug for it returns an 
erroneous matrix. The cause of this weird occurrence remains unknown, however we 
shall attempt to effectively overcome this mishap by applying an easy trick. 
i= 0. rows(new_datamatrix02) — | 
ro cols(new_datamatrix02} 


cor_norm_diff = submatrix{ norm_diff,, 0,rows(norm_diff, — 1,0,r. - 1 


This last matrix named "corr_norm_diff" is the one we intended to find. In each row of 
unequal element number, it contains a different random variable as function of our 
"devised" groups. Let us take a closer look in one of them, say the 5th. 

| \\O? 
ext (tone mune 


} = -0.156 


mean(test_O p 


oe 


hi 


test_ 1 <= (oe Nore: 
mean(test_l = -0.059 
Detining the number of the bins: 
Nnoints = 30 

j= 0.. rows(corr_norm_diff) - | 
Bato ‘= max corr_norm_dit 

oS = mir{ corr_norm_diff 


upper_max := man upper ) 
lower min:= minclower) 


UpPehlindx = ower niin 


SC = 
N 
points 
ven_intervals roa lower_min+ step-1 
hist ( Npoints »COrm_norm_diff fl 
L := 


J cols(corr_norm_ditf ) 
J 
step = 0.071 
buildL(L) := |Lnew — Ly! 
(Omee Pa pows (le). 


Ae¢L 
q 


Lnew — stack eacas A) 


Lodt = DUNGLANS) 
rows( Lay) = 12 
cols(Lygr} = 30 


i.) 


> boa, =! 


— 


=) 


13S 





Lodt 


ave _pdt = ee | 


 avg_pdt = | 


mean(avge_pdf ) = 0.033 
Defining the Gaussian pdf: 
mean_pdf ‘= mean(avg_pdft ) 
mean_gauss := Q-mean_pdt 


Taverod? > 
var_pdf := val) ——— | 





step 


/ 


var_gauss := 0.65 var_pdt 


Sidev anal =.= y var_pdt 
stdev_pdf = 0.42 
var_gauss = 0.115 
= /gen_intervals,—mean_gauss} ~ 


2 var_gauss 


las) 


URTUSS te 
i 1] 
2-7 var_gauss 





Ibi, 


U | 


0 0S 


0.06 


0.04 


0 02 





0 

| =e: 0 0.5 | 15 
*€ average pdt 

==" gaussian pdt 


, 


UBoun “= 20-log( I + gen_intervals ‘| 
UBioss °= —dBeoain 


0. | 


0.08 


0.06 
uv pal, 
UUs) SLED 
r= 0.04 


002 
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1] 
is 1 
m=O 


0.8 
0.6 


0.4 





; a 
= 1G =) 0 
OK cdf 


S27 Se oauisstan- cal, 


LAN 
S 
SAI 
Io 
(=>) 
tt 
ta 


End of typical room example (random variable calculation) 


MathCAD file used in the room with furniture example (No PEC backing) 
Reading the setting the data: 


cls >.99.10° 
f-=2)000 10° 
—— = 

, 
A = 0.299 


For the PEC backing case we simply change the data files. 
BE = READER IN BSC resultsOS.tx1. } 

SBR := READPRN SBE _results.txt” ) 

EFree = READPRN"BSC_results_freeQ2_1.txt" ) 


PLsBR dB = SBR’ Y 


0, rows(PLegr ap) — | 
27 = SBR. 


3 
Eig ke 


Eydp (= EE 


ee 


ion 




















E,dp = EE 
at 
EL ap 
20) 
Ek, = 10 
a 
E. ap 
20 
Ey =O 
—_——>) 
E7dB 
20 
Ey and 
G)) 
Ey treedB '= EGree 
es. 
EY ireedB «= EEfrec 
EytreedB = EEree 7 
; 
EE frcede 
20 
Extree = 10 
oe 
EL freedB 
- 20 
Ex tree = 10 
——— 
Seeds 
20 
Evtree = 10 





2 2 2) 
Frot_tree “= Exfree + Eyfree + Ezfree 


2 2 2 
Bor = Ex + By” +E 


iVe= 0 vrows(z) — | 


[ ‘e | ae ] 


Plier = 1Olvg 
at] 
n 


Calculating the excess path loss and the NEC-BSC predicted path loss 


——— 


Etot 


Excessyp := 20-log 
Exot_tree 


PlasC = PFliree Px cestap 


x 


rove F501 Oo. 
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Comparison: SBR and NEC-BSC w/o PEC 


~40 


Path Loss (dB) 


—70 





distance (m) 
ae Gy Se sates 
ee SS Bi alneoss 


PLasc 
; 1O 
Plaew = 10 
Nn 
PLspr uB ; 
10 
pilspr :=10 
m 


——9 
i 
mean_square_BSC := net plan? | 
-9 
meanosatare BSG = 6.003 10 
| | 
mean_square_SBR := mean my 


mean_square_SBR = 1.866x 10 : 
mean_BSC:= mean(plnew) 


mean_BSC = 4.26x 10. : 


mean_SBR := mean(plspp) 


5 


Meus BR = 71/37 x 10- - 

SBR_power := meun_square_SBR 
SBR_power = |.866~x 10 ; 
BSC pO Wer = meanosquare. BSC 


Loe 


- 9 
BSC_power = 6.003 10 
sel_ :=m20 

m 
BSCcomp_ = Plnew cel 
rows| BSCcomp) = 70 


diff fGerEC y= BSCcomp - pISBR 
1) in 1) 


WRITEPRN "no_PEC.dat” ) ait no_PEG 


dit_no PEG]. 





2Z m 


End of room with furniture example (no PEC backing) 


MathCAD file used in the room with furniture example (comparison) 
Compuring the no-PEC and PEC results 

with PEC = READPRN “with_PEC.dat” ) 

without PEC *=READPRN no _PEC dat) 

m:=0.. rows( with _PEC) - | 

fee Se in 


Comparison-with & without PEC backing 


peer Ee aAeg P= 
a0 Op Petdd. 


Bes 8 & A 4 
<p: ° oe. a 
4 


FN oN eo ein te ANNAN Bas 
4 71 ae sn Ne a9 <t Ca eee 4 
& 


difference (linear) 





distance (m) 
IK with PEC 
+--+ without PEC 


Cal Within Go) = 2 kOs x 10. : 


~ 


mean( with _PEC) = 1.239x 10 . 


Vah wtnOulue ee Ce= 1.05 5>< 10. : 
-5 
meant without_PEC) = -3.149x 10 ~ 
[ 2 -§ 
mean| with_PEC } = 2.119x 10 
er 
( a -§ 
mean| without_PEC™ } = 1.154 10 
7 — 


S 

| meant with PEC | 

lOlog ae a 
{ \ 


2 
\ mean| without_PEC 


End of room with furniture example (comparison) 


= 2.638 





MathCAD file used for the calculation of the composite walls constitutive 
parameters 


Computing the relative permitivity for composite walls 


In the input NEC-BSC file we come across basically two composite walls being located 
close to antenna’s room. The first one ts not a real wal] but simply a meta] mesh which we 
model as an “1maginary wall.” The second one consists of approximately 80% plaster and 
20% metalic studs. 


G5 


The formula dernved last ume for the mesh was the following: 
1 wa af e r 
hes a Er 


Z 


| PS fice _ metal_area — a 

eee Jer ~  wall_area A 
(equauon 1) 
where “epsilon subscript r” is the relative permitivity of the “imaginary maternal” 
composing the mesh. Bear in mind that if A is the total surface of the mesh and a the total 
surface of the metal structure then a/A equals 0.17355. The above equation 1 does have a 
real solution. 
Composite material case I 

We consider the following approach. An EM wave with power Pin impinges on a 
wall consisting of plaster and metalic studs in proportion 80 and 20% respectively. The 
transmission coefficient of this wall 1s Tcomp. We then let Tplaster being the 
transmission coefficient of a plaster wall with surface 0.80 times the wall’s surface. It 
must be true that: 


2 
Ptransmitted= (1 — 0.2) Pin 1 a ( iprasten | 
(equation 2) 


E 


Ptransmitted = Pin 1 7 eran fal 

(equauion 3) 

REFLECTED POWER IS PROPORTIONAL TO THE MAGNITUDE SQUARED 
OF THE REFLECTION COEFFICIENT! 

assuming the transmission due to the metalic portion of the composite wall being 
negligible and therefore contributing negatively to the resulting transmitted power. 


For normal incidence, we have: 


J -./ €r_plaster 


Pilster = = 
1+. €r_plaster 


€r_ comp 


ee 
[ comp = 
1+ €r comp 


Er_plaster = 2-3 
(equauions 4) 
Thus, from equations 2,3, and 4 we obtain: 


<3 
1-V€r comp | I -./€r plaster 
Ye 4 ee EV p= 108 ee eee 
eee | i+ €r_ plaster | 
(equalion 5) 


trom which it follows that: 


i) 
ny 





r_ comp 


} 
4 ~ 
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mS = sco xt l= 0 
(equation 6) 
where x equals (epsilon_comp)(1/2). Therefore: 
x1 = 2.8717 
or 
x2 (= 0.3483 
(equation 7) 
And finally: 
5 
€r_comp1 ‘= *1 
or 
2 
Er_comp2 ‘= *2 
€r compl = 8.247 
Er comp2 = 9.121 
(equauons 8) 
Composite material case IT 
In this case we consider the metalic mesh being "replaced" by an imaginary wall 
consisting of air and metal in proportions 82.6447 and 17.3553%, respectively. We 
assume that only the “air-part"” of this imaginary wall contributes to transmission and only 
the "metalic-part"” reflects back EM energy. Therefore, if Twall and Rwall are the 
transmission and reflection coefficients, respectively, of the "imagnary wall” then (no 
absorption occuring): 


Zs 
Prransmitted = ( 1 a O5173533) Pin| 1 a (lve) | 
(equation 9) 
2 
Ptransmitted= Pin iF 7 (oes 


(equation 10) 
where Gamma_mesh equals: 


rf 1—./€r mesh 
inci —= == 
hee / © r_mesh 


und 
Dair = 0 
(equation 11) 


From equations 9,10 and 11 it follows that: 
9 


I / © r_mesh 
1+,4/€r_mesh 
(equation 12) 


From which we obtain: 


0.826447 =| 1 - 


}o 
O? 
~] 


e —-2.84-x+1= 0 
(equation 13) 
Thetetore: 

X] (= 2.4281 

rx» '= 0.4118 
(equation |) 
And finally: 


fr meshl -= *1 


Er _mesh2 := X2 

Er meshl = 5.896 

or 

Er _mesh2 = 9.17 

(equation 15) 

End of MathCAD file used for the calculation of the composite walls constitutive 
parameters 


MathCAD file used in the comparison with measurements in Bullard 2nd floor 
Comparison between measured and predicted results in Bullard 


Reading the measured data 

hall-= READPRNG puller |) 

rows (hall) = 29 

Colt hall) = 24 

Cropping in order to rule out the "ambiguous" regions. 
cropped_hall := submatrix (hall, 0,28, 4, 23) 

rows (cropped_hall) = 29 

cols(cropped_hall) = 20 

Reading the NEC-BSC predicted data 

Hredivicd .= READER NG prediciedixt | 

rows (predicted) = 696 

vols¢ predicted) = 9 

predicted *’ | 


a | cor 


> 
(5) 
predicted *~’ 


R210 7 
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— 


predicted 7) 





rOWSsS Bion = 696 


Converting the values to dBm 
Gdipole = 1.64 


Gmonopole ‘= Gdipole 
Gmonopole = ioe 


= 0.122364 
ma 
A ‘Gmonopole 
Aettective “= 
4-T 
Aettecuve = 1.954x 10 
eee 
2 
Eel) 
density - 2.377 
ee 


Preceived ‘= Pdensity ‘Affective 

Creating the NEC-BSC data matrix in the same fashion with the measured data matrix. 
N.x is the # points in the inner loop 

\ y is the # points in the outer loop 

ety ao 


EN ee, oe 





y= 
Nx 

Nad 

ed) ae | 

room * := Submiathi «i eeses ced . Slarly , Stopg .0, 0) 
—— 

Preceived_dBW — KOloe lreneived 

room. dBW’ = submatrix(Preceived dBW , Staltg , SLOPg 0, 0) 


rOWws [Ree cdnae = 696 
rows (room_dBW) = 29 
cols(room_dBW) = 24 

mux ToOmece WwW) = =1451 
min¢(room_dBW) = —251.093 


Bohs, 


maux_power_dB := max(room_dBW) 
dyn_range := 30 

power_limit := max_power_dB — dyn_range 
power_limit = —44.51 


p :=0.. rows (room_dBW) — 1 


room_power_dBW), ‘= if(room_dBWp g > power_limit, room_dBWp g, power_limt) 


( - 
room_lin) := submatrix( Preceived > Starlg, Slopg . 0, 0) 
rows (room_lin) = 29 


colsiroom_ lin) = 24 





++_dBW 


Converting to dBm and "cropping" appropriately. 
oe 

Ptranspose_dBm := room_power_dBW + 30 

Preceived_dBm?= Mig es ee Sa 

min( Ptranspose_dBm) = ~—14.5] 

max(Ptranspose_dBm) = 15.49 

rows (Preceived_dBm) = 24 


Precel\ ed dBm) = 529 


}- = 


rows (Ptiranspose_dBm) = 29 

cols(Ptranspose_dBm) = 24 

cropped_Prec_dBm := submatrix(Preceived_dBm, 0, rows (Preceived_dBm) — 1,4,cols(Preceived_dBm) — 1) 
cropped_Ptrans_dBm := submatrix( Ptranspose_dBm, 0, rows (Ptranspose_dBm) — 1,4, cols(Ptranspose_dBm) — 
rows (cropped_Prec_dBm) = 24 

cols(cropped_Prec_dBm) = 25 

rows (cropped_hall) = 29 

cols(cropped_hall) = 20 

rows (cropped_Ptrans_dBm) = 29 

cols(cropped_Ptrans_dBm) = 20 


Converting to linear in mWatts 
a 
cropped_hall 


linear nall:= (0 ~ 


max_linear_hall:= max(linear_hal} 
max_linear_hall= 3.162x 10° 


——_—_—_—_—_—_—— 


cropped_Ptrans_dBm 


hnéear Pree -— 0 ” 


max linéaperrces—= max (linear Prec) 
max_lin@amee fee ——9 451 

Scaling accordindingly 
linear_Precjg. 14 = 0.072 


| 

max oimeat halls cue. ee 

ol w ss  ineal ar rer 
max_limear_Prec 


scaled_BSCig. 14 = 24.251 


Converting the scaled BSC values to dBm again 
SEE 


BoGsdbin-— 10 loc(scaled. bSO)) 


Calculating the difference between BSC and measured values 
: faa generis eae 


difference_dB := (cropped_hall - BSC_dBm 

difference := (|linear_hall- scaled_ BSC] ) 
meian_dift_linear:= mean( difference) 

meun_ditf_linear= 245.294 

in mWatts 

Reversing the matrices to improve visualization 
qq := 0.. cols(cropped_hall) — } 


5 H)-1- 
ee er = cropped_hall°o's(eropped_ha ) qq) 


yeuhecel ane eeu = BSC_dB ;\c0ls(BSC_dBm-1-gq/ 


jug eat 


BSC_dBm 
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yi elitial 
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measured_map 





predicted_mup 


End of MathCAD file used in the comparison with measurements in Bullard 2 
floor 


MathCAD file used in the comparison with free space (Bullard 2nd floor) 


Piel Or ihe Boek eculic 

Read the Data for antenna at location |. For the rest of locations we simply change the 
duta tiles | 

room_data := READPRN("CH6_onginal_position.txt 

free_data:= READPRN("BSC_Bullard_free_space.tx)’ 

M := rows(room_data) 


GOA 10° 


cols(room_ data) 


M 
hae 


Lis 


N=Y 


phi := feeipwdeg as 
eae 
E = \ i a 
ee! 
gh aaa 
? 
Fo = 10 
a 
Pe eae 
> 
Eg := 10 ey 
free datas” 
*) 
e; =| 10 a 
free_data™” 
a 
eg = 10 ee 
a 
A 
Co = 10 A 
ee a 
2 2 ws 
Ew = y (Er) + (Eo) + (Eo) 
ee 
wi 2 Zz 
erot = yf (er) + (e6)° + (£0) 
n= 0 Shi =) 


N.x 18 the # points tn the outer loop 


row s(een = 2.624 x 10" 





de-N 4 


Slarlg <= q-Nx 


174 


StOPg := Slarig t+ Nx — | 


on = submatrix Ej, . $larlg , StOPg, 0, 0) 


-SSSraS 
Erot_ dB = 20-logq{ Exot) 
Ge 


room_dB := submatrix Exot_ dB. Statlg, Stopg 0,0) 

maxioomecs)) = 42°59 

min room_dB) = —195.229 

iid ile lms —inax Foon cb) 

dyn_range:= 30 

field_limit= max_field_dB- dyn_range 

p := 0.. rows(room_dB) — 1 

room_fiel@ |g := if(room_dBp 4 > field_limiroom_dBp , q,field_lim)t 


nen = submatrix Exot , Starlg, Stopg 0,0) 


free? = submatrix eto, : Starlg ; StOPg : Ue 0) 


—— 
Ctot_dB “= 20-log{etor] 


free_dB = submatrix etor_dB , Slarlg , StOpg 0, 0) 


max freee) Se2 | 

minmuiee db) =—s.7 / 

max_field_dB= max(free_dB) 

dyn_range:= 30 

field_limit= max_field_dB- dyn_range 

itelealinntity oo | 

p := Q.. rows(free_dB) — 1 

free_field q := if(free_dBy g > field_limitfree_dB)  q,field_lim} 

free lin? = submatrix err , Stat , StOPg , 0 , 0) 

diff_free_room= free_lin- room_lir 

diff_free_room_dB = if( diff_free_roorg ,g > 0, 20-log{ diff_free_roormp, g) 0) 
Mil vil = aici etree TroOmmdlp 

min_val= —92.194 

diff_fr_dB q := if(diff_free_room_dB q = 0,min_valdiff_free_room_dB, q) 
rows(room_field = 162 

cols(room_field = 162 


Wes 


ir = 0. rowsGoon eda 


.) ] panes ee |. 
nversed: romance = room_fiela 7S (room_tield) ) I-11) 





a el" 


a le a 
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inversed_room_field 
End of MathCAD file used in the comparison with free space (Bullard one floor) 


MathCAD file used in calculating the excess gain/loss random variables 


In this program we are aiming at calculating all the random variables of the excess 
gain/loss of the room over free space field strength (Etot and Etot_free respectively). 
We start by reading the associated data files generated by NEC_BSC. 
NEC_BSC_data:= READPRN "CH6_original_posiuon.txt” ) 

M := rows( NEC_BSC_data) 


4 
ME = 22624 1) 


= ColsNec BSC daia) 
eee 


(Q) 
‘coord = NEC_BSC_data 


} 


P= NECeBS ead. 


\ 
« COO! 


Satake (9) 
Zcoord : NEC _BSC_ data 


Antenna coordinates for Location |. For the rest of locations we simply change the 
coordinates: 


t =| 
= 
oO’ 


Kant c= 20 

Yant °= 8.0 

Zont i= | 

Distance from the antenna: 
ee 


9) ) 7 
Pp i= beard = Xant) + eae a venn + fear = Zan 
free_space_data := READPRN"BSC_Bullard_free_space.txt" ) 
MM := rows(free_space_data ) 
MM =2.624x 10" 
NN := cols(free_space_data ) 





NN = 9 
, 
. (3) 
NEC _BSC_ data 
20 
Ey 10 
es 
NEC_BSC_data™ 
20 
ee as NG 
ne 
NEC_BSC_data’” 
20 
0 
es 2 2 
Erot =f (Ex) + (Ey) + (E,) 
4 
rows(Eior ) = 2.624 10 
cel spneemnn. 
20 
Ex free “=| 10 
ae 
fee soacenae 
20 
Ey tree = 10 
ee 


re 


20 
E, free = 10 


SSS ey 
y 


2 2 2 
ELot_tree = Meer! a Tey hee) i ice! 
Since the NEC-BSC dynamic range is only 30dB, we will have to account for this fact 
by “filtering” out all values which are out of that range. This is done in the following 
fashion: 









——_—_— 


Erot_ dB = (20-log(Eror )) 


her 


dvn_range := 30 
=U. rowslEn ap) ! 
lowest allowed Etot_value *= man Eo dB) * dyn_range 
lowest_allowed_Etot_value = 12.59 
8 a 
Frot_free_dB ‘= (20-log(Exor_tree)) 
lowest_allowed_Etot_fr_val := man Erot free dB) — dyn_range 


lowest_allowed_Etot@ipavale= 702! 
Eotal dB. = if (Ero dB F < lowest_allowed_Etot_value , lowest_allowed_Etot_value ,E,or ap | 
aoe ee ~ 


Ee Ainlaive: 4k {Exot tree dB. < lowest_allowed_Etot_fr_val . lowest_allowed_Etot_fr_val ,Exor free aB | 
iT aoe | ee 
Converting to linear: 
ee 
CLotul_dB 
20 
Exotal ‘= 10 
— — —-- 
E otal_free_dB 
20 
Fotal_tree = 10 


Let us now create the two datamatrices, namely the “datamatrx" which consists of the 
rooms field strength values for every observed point and the "freespacedataO1"” which 
includes the corresponding data for the tree space case. 

datamatria := augment(p doen) 

treespacedataQl := augment(p Bera isitce) 

At this point, owing to the non-uniform layout of the distance (from the antenna) points in 
the _ array, we will have to distinguish them into several different groups using the 
following approach: 

lower= floorl min p )) 

upper ‘= ceil maxp)) 

Nlaximum and minimum distance for this antenna location: 


upper = 49 
lower = 0 
Groups separation for this calculation (in cm): 


separauion := 50 
"Step factor” for this calculation: 


100 | 
tucor =cel ———_— 
_ separauion 
ee (ole 


Number of distance bins for this calculation: 


bins := upper -factor 
bins = 98 


BBS: 


7 upper - lower 
bins 

h=O5 

i = 02 bins 


K ‘= lower + n-h 


[t should be noted that the number of bins may be chosen in a fashion which best suits our 
desirable “distance steps’. For instance, should we choose to have groups points 
separated by 50cm. where the maximum distance from the antenna 1s 49m, our intuitive 
selection should be 98. Unfortunately, for reasons that could not be explained at the trme 
of thesis write-up, the “build” function does not yield correct results for some values (e.g. 
for separauion equal to 25cm). 

We next introduce a new function, which shall be utilized, several times in our 
subsequent calculations. As input, it expects a two-column matrix and an array. It looks 
for those values of the first column (1.e., in our problem, the distance array J 0 Owhich 
fall within our made-up groups and then creates a new matrix that contains all those 
second column's values whose respective first column value satisfy our requirement. It 
finally returns this new matrix, which, it should be underlined, normally has a different 
number of elements in each row. 


(0) 


build(M.binvector) := }A — M 


Be vi 

k < binvector 

for le 0.. rows(k) -— 2 
D, < 0 

LOE whe Or-Taws(k) == 2 
for j;€ O..rows(M) — ] 


Dee augment(B..D | ieee See 
| : Jol J I+ | 


tor 1€ Q..rows(D) - | 
D = submatrix{ D,,0, rows(D.) 2 1,0, cols(D,] 7 2) if D. #0 


D 
It is high time we used this helping function. 
datamatrix02:= build(datamatrix, k) 
treespacedataO2 := build(freespacedataQ1! , k) 
The following new function 1s, seemingly, of high importance since it 1s the one which 
calculates the "normalized difference" random variables. It receives two matrices of the 
same dimensions and returns a new matrix, which in each row contains the values of the 
random variable pertinent to that row’s "distance group." 


ves 


difference( X.Y) := | for ke 0.. rows(X) - 1 
Cc, —0 
for j}€ 0.. rows(X) - 1 
ne cols( X 
Ac = 


ee 





i1— 0 
for i€ O.n-] 


A 
| diff — 20log 





ee 
ne cols(X,) 





! Cie stack(C.. diff] 
: J J 
c= (c) 
ke rows(C 
le cols(C.) 


| i) oe submatrix{ C. ,0, k =e eee 1} 
J J 


C 
The "room gain” is defined here in dB 


We now get rid of the redundant first row zeros and then calculate the normalized 

difference variables. 

new _datamatrx02 := submatrix(datamatrx02 1, rows(datamatnix02 — 2,0, cols(datamatrix02 — 1) 
new_treespacedata02 := submatrix(freespacedataQ2 , | ,rows(freespacedata02 ) — 2,0, cols(freespacedataQ2 ) - 1) 
norm_ditt := ditference(new_datamatrix02, new_freespacedataQ2 ) 

Surprisingly and unexpectedly Mathcad appears to suffer form a bug for it returns an 
erroneous matrix. The cause of this weird occurrence remains unknown, however we 

shall attempt to effecuvely overcome this mishap by applying an easy trick. 

1= 0.. rows(new_datamatnxQ2) - | 


= cols/ new_datamatrix02] 
| 


con_norm_ditt = submatrix{ norm_dift,, 0, rows{norm_ditf,) - 1,0,r, — 1 


This last matrix named "corr_norm_diff" is the one we intended to find. In each row of 
unequal element number, it contains a different random variable as function of our 
“devised” groups. Let us take a closer look in one of them, say the 5th. 

\4O? 


(est: = (Ganache 


mean| test_0 5] Sod 
7\ 

(si = ermaormeditt 
meant test_1 0) =i 006 
Defining the number of the bins: 
Npoints ‘= 500 
ie OO comunoimediti) a) 
upper := ma ‘corr_norm_ditt | 

pper 
ic = min(cor_norm_dift 
upper_max := max upper ) 
lower_min:= min(lower) 
upper_max — lower_min 
Slep ‘= oa ae 

Npoints 

upper_max = 25.861 
lower _min = —14.794 
siep = 0.051 
ven_intervals — lower_min+ step-1 


hist ( gen_intervals , com_norm_diff 


C= - 
J cols (corr_norm_ditf 


buildL(L) := |Inew — Ly. 


for qé l..rows(L) - | 
AeL y 
q 
Il © |Enew <— stack(Lnews A) 
Lop := buildL(L) 


rows(Lpgf ) = 96 
cols(Lyar) = 499 

Npdfs = cols(nar} = | 
(Ge Ue Npdic = 


| (aq 
avg pdt a = mean| Log a | 


Save pat = | 


N pdis7! 
mean_rv := > gen_intervals -avg_pdt 
qq qq 


qq =0 
mean_rv = 7.165 


Leo 


step = 0.081 


‘\ | 


Pus” 


o 
Mean_sq_ rv = » (gen_intervals || "-avg_pdt 


yy = 0 
mean_sq_rv = 75.007 


P 
var_rv := mean_sq_rv — mean_rv_ 


“olen = Moved! 
mean_rv = 7.165 
Defining Gaussian pdf: 


4 


- (ge n_intervals;—-mean_rv) 
| 2: var_rv 
2aUuss = 
2-1: Var_fv 


room_gaingp ‘= gen_intervals 


Antenna at location | 


0.08 


0.00 


pdf 


0.04 


Owe 





ss 10 he) 0 5 10 [5 20 25 30 
Room gain (dB) 
<K average pdf 


== ° gaussian 
qq 
hae — yy avg_pdt 
qq m 
m=0 
NI 
cdf gauss, = step | »; Qauss m | 
| im =o j 


Antenna at location | 
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End of MathCAD file used in calculating the excess gain/loss random variables 


End of all MathCAD files 
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